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I HlSTORIQUE 

Au debut de l'informatique, le fait de programmer se revelait fastidieux et peu fiable. Puis apparut 
une codification directe des instructions machines, connue sous le nom de code assembleur. Dans les 
annees 50, il devint evident qu'une nouvelle etape dans les techniques de programmation devait etre 
franchie. 

C'est pourquoi, l'annee 1954 voit naitre le projet du premier langage symbolique par John Backus 
de la societe IBM ; Mathematical Formula Translating System dit FORTRAN. Ce langage permet pour 
la premiere fois a l'utilisateur de se concentrer sur le probleme a resoudre et non sur les instructions de 
la machine. 

1956 - Premier manuel de reference qui definit le Fortran I. Les noms de variables ont jusqu'a six 
caracteres et l'instruction FORMAT est introduite. 

1957 - Le Fortran II fait son apparition ainsi que les premiers compilateurs commerciaux. Les 
sous-programmes et fonctions sont acceptes et deviennent compilables independamment les uns des 
autres. 

1958 - A la fin de l'annee, le Fortran III est disponible, mais reste interne a la compagnie IBM. 
Certaines des nouveautes implantees sont le type logique, les parametres de procedures. 

1962 - Le Fortran IV est ne et restera le langage informatique des scientifiques pendant seize ans. 
II apporte les declarations de type explicites et l'instruction d'initialisation DATA. Entre temps, les 
autres constructeurs developpent leur compilateur ; il faut bientot etablir une norme. 

1966 - Creation de la norme ANSI (American National Standards Institutes), le Fortran IV est 
rebaptise Fortran 66. 

1978 - La norme Fortran 77 (ou Fortran V) modernise le langage, car c'est la fin des cartes 
perforees. On lui doit le type caractere, de meilleures entrees / sorties, un debut de programmation 
structured et des fonctions predefinies portant un nom generique. 

Dans le milieu des annees 80, la reconversion a Fortran 77 etait bien engagee alors que Fortran 66 
perdait considerablement du terrain. 

Apres 30 ans d'existence, Fortran n'est plus le seul langage de programmation disponible, mais il 
est reste predominant dans le domaine du calcul numerique a l'usage des chercheurs, des techniciens et 
des ingenieurs. Dans le but de maintenir le langage correctement a jour, le Comite technique X3J3 
accredite par 1ANSI a elabore une nouvelle norme. 

1991 - La norme internationale ISO/ ANSI Fortran 90 apparait dotee d'innovations changeant 
profondement le langage. Le Fortran devient modulaire, il permet un calcul vectoriel et un controle de 
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la precision numerique. De plus il devient plus fiable grace aux blocs interfaces et s'oriente vers une 
programmation objet. 

1995 - Les caracteristiques depreciees ou obsoletes forment le Fortran 95 et il admet une structure 
vectorielle supplemental (FOR ALL) et quelques autres ameliorations. 

2002 - Fin des travaux pour la norme Fortran 2000, publication prevue pour fin 2004 pour la 
norme ISO. 

Fortran 2005 : Correction de Fortran 2000 et probablement derniere norme Fortran. 
Plus tard, va-t-on vers un langage F, FOO, F++ ? 
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II 



Structure d'un programme et syntaxe 



II.A 



Unites de programme et procedures 



De facon tres generale, la taille des programmes que Ton est amene a developper aujourd'hui 
devient considerable. Les ordinateurs sont de plus en plus puissants, ils sont alors de plus en plus 
sollicites dans la vie courante. Et l'exigence de l'utilisateur n'a plus de limite. Quant aux chercheurs et a 
leurs desirs, rien ne saurait les decourager, pas meme des dizaines de milliers de lignes de calculs. Or, 
des que Ton envisage un travail un tant soit peu ambitieux, on a de plus en plus sou vent besoin de 
competences multiples. C'est-a-dire que les programmes sont developpes par plusieurs personnes, qui 
hormis le langage informatique commun ne parlent pas forcement le meme jargon technique. 

C'est pourquoi, nous avons besoin d'un langage structure et d'un langage modulaire. La structure 
permet une bonne lisibilite d'un programme, done une bonne comprehension et une maitrise du codage. 
C'est-a-dire que Ton resout effectivement le probleme donne et pas un ersatz qui fonctionne dans 97 % 
des cas, voire moins ! La modularity permet de repondre aux besoins du Genie Logiciel. Les modules 
fonctionnent quasiment comme des boites noires, ce qui permet d'incorporer dans un programme des 
parties independantes les unes des autres et ne necessitant pas forcement aux autres developpeurs de 
prendre connaissance de leur contenu. II suffit de connaitre le nom d'appel des procedures dont on a 
besoin. 

La notion d'objet est plus forte encore, et on les manipule alors completement comme des boites 
noires. L'objet realise lui-meme l'encapsulation de ses donnees et empeche tout acces depuis 
l'exterieur. 

En recapitulant, Fortran 90 facilite la manipulation des structures (objet de type derive), en offrant 
les possibilites suivantes. 

- II permet de definir des fonctions de manipulation de ces structures : 



- II permet d'encapsuler le type et les fonctions operant sur les objets de ce type dans un module 
separe. Ce qui permettra d'y acceder facilement (USE <module>) en diminuant les sources 
d'erreurs, dans toutes les unites de programme en ayant besoin ; 

- II permet de rendre privee les composantes de la structure. Si la definition du type se fait dans 
un module separe (encapsulation), seules les fonctions ou procedures internes au module ont 
acces aux composantes privees (PRIVATE) et peuvent les modifier (securite). 

Sans etre un vrai langage oriente objet, Fortran 90 fournit ainsi des extensions objet bien utiles 
pour la structuration du programme et sa fiabilite. 



- surcharge des operateurs +,-,*, 

- de nouveaux operateurs en les associant aux fonctions correspondantes, 

- autres fonctions ou procedures agissant sur ces structures ; 
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Un programme est compose d' unites de programme, qui sont : 

- un programme principal, contenant une procedure interne ; 

- des procedures externes ; 

- des modules. 

Les procedures externes peuvent etre ecrites dans des modules, elles deviennent internes au 
module. Un module peut contenir des procedures, des declarations, des affectations de variables, des 
types derives et des blocs interfaces. 

Les mots cles sont : 

PROGRAM, END PROGRAM ; 
MODULE, END MODULE ; 
CONTAINS ; 

SUBROUTINE, END SUBROUTINE ; 
FUNCTION, END FUNCTION. 

II. A. 1 Le programme principal 

Le programme principal est une unite de compilation, mais il a la particularite d'etre la seule qui 
soit executable. Done il est la seule qui soit totalement obligatoire. Le programme principal prend 
toujours la forme suivante : 

PROGRAM nom_du_programme 

[instructions de specification] 

[instructions executables] 
CONTAINS 

[procedures internes] 
END PROGRAM nom_du_programme 

L'instruction PROGRAM est facultative, mais conseillee. L'instruction END qui genere l'arret du 
programme est obligatoire. Les instructions de specification sont des instructions de declaration, des 
definitions de type, des blocs INTERFACE. L'instruction CONTAINS indique la presence d'une ou 
plusieurs procedures. 
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Exemple II. 1 

PROGRAM ETOILE 
USE DONNEES 
IMPLICIT NONE 
REAL :: rl, ul, Tl 

rl = l.e-3 ; ul - 45000. ; Tl = 55635. 
rl = G * ul * Tl 
PRINT*, "rl vaut ", rl 

END PROGRAM ETOILE 

II.A.2 Les procedures externes 

Une procedure externe est appelee depuis une autre unite de programme. Elle a une forme 
similaire a celle du programme principal, tant du point de vue de sa syntaxe que de son execution, elle 
prend place On en distingue deux sortes : les sous-programmes et les fonctions. De plus, on distingue 
des procedures globales, contenues dans un module et utilisables par l'ensemble des autres unites de 
programme, des procedures internes, contenues dans une procedure globale et utilisables uniquement 
par leur note. 

a) Les sous-programmes 

Les sous-programmes, definis par le mot-cle SUBROUTINE, definissent une action au moyen de 
toutes les instructions du langage, dont l'appel prend lui-meme la forme et l'apparence d'une instruction 
du langage. La syntaxe generale d'un sous-programme est la suivante : 

SUBROUTINE nom_du_sous -programme [(arguments)] 

[instructions de specification] 

[instructions executables] 
[CONTAINS 

[procedures internes]] 
END SUBROUTINE [ nom_du_sous -programme ] 

Les sous-programmes sont appelees, depuis une autre procedure (ou programme principal), avec 
l'instruction CALL nom_du_sous-programme [(arguments)]. 

b) Les fonctions 

Les fonctions, definis par le mot-cle FUNCTION, calculent une valeur au moyen de toutes les 
instructions du langage, dont l'appel prend place au sein d'une expression. La syntaxe generale d'une 
fonction est la suivante : 



10 



Ecole doctorate d'Astronomie et d'Astrophysique d'tle-de-France - Cours de Fortran 90 1 95 - Claire Michaut 

[TYPE] FUNCTION nom_de_la_f ontion (arguments) 

[instructions de specification] 

[instructions executables] 
[CONTAINS 

[procedures internes]] 
END FUNCTION [ nom_de_la_f onction ] 

Les fonctions sont appelees directement a l'aide de leur nom propre et de leurs arguments, depuis 
une autre procedure (ou programme principal). On peut declarer leur type avant l'instruction 
FUNCTION. 

Exemple YY.2 

SUBROUTINE LOCATE (rr, TT, m, n, hmn) 
IMPLICIT NONE 

REAL, INTENT (IN) :: rr, TT 
INTEGER, INTENT ( OUT ) : : m, n 

REAL, OPTIONAL, DIMENSION ( 4 ) , INTENT ( OUT ) :: hmn 
REAL:: hm, hn, hml , hnl, Vol 

m = 1 ; n = 1 
DO 

m = m + 1 
n = n + 1 
END DO 

IF ( PRESENT ( hmn ) ) THEN 

! Calcul des poids pour 1 ' interpolation des donnees. 

hm = TT - TeV(m) ; hml = TeV(m+l) - TT 

hn = rr - rho(n) ; hnl = rho(n+l) - rr 

Vol = (hm + hml) * (hn + hnl) 

hmn(l) = hm * hn / Vol 

hmn(2) = hml * hn / Vol 

hmn(3) = hm * hnl / Vol 

hmn(4) = hml * hnl / Vol 
END IF 

END SUBROUTINE LOCATE 
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Exemple II.3 

REAL FUNCTION ZMOY (rr, TT) 
IMPLICIT NONE 

REAL, INTENT (IN) :: rr, TT 
CALL LOCATE (rr, TT, m, n, hmn) 

zmoy = zbar(m,n) * hmn(4) + zbar (m+1 , n ) * hmn(3) & 

+ zbar (m,n+l ) * hmn(2) + zbar (m+1 , n+1 ) * hmn(l) 
zmoy = 1 . + zmoy 

END FUNCTION ZMOY 

c) Les procedures recursives 

Dans certains cas particuliers, on a besoin qu'une procedure puisse s'appeler elle-meme. 
D'ordinaire, cette structure de programmation est impossible. C'est pourquoi, il faut tout d'abord avertir 
le compilateur de ce mode de fonctionnement : on utilise l'attribut RECURSIVE devant le nom de la 
procedure, qu'elle soit un sous-programme ou bien une fonction. Attention toutefois au temps de calcul, 
qui peut devenir prohibitif, lors de plusieurs appels recursifs successifs ! 

Syntaxe generate : 

recursive subroutine nom_proc (list_arg) 

recursive function nom_proc (list_arg) result (nom_res) 

Remarque : dans le cas d'une fonction, on est oblige d'utiliser un nom de resultat different du nom 
de la fonction, defini a l'aide du mot cle RESULT. 

Les fonctions recursives peuvent etre, par exemple, tres utiles pour traiter des suites 
mathematiques, qui par definition, sont construites en reference au rang inferieur. Ainsi pour calculer la 
suite arithmetique : 

U\ = \\ u-2 = 3 ; u n = 3u n _\ - u n _2 
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Exemple II.4 

RECURSIVE FUNCTION SUITE (n) RESULT (u) 
INTEGER, INTENT (IN) :: n 
INTEGER : : u 
SELECT CASE (n) 
CASE (1) 

u = 1 
CASE (2) 

u - 3 
CASE DEFAULT 

u = 3 * suite (n-1) - suite (n-2) 
END SELECT 
END FUNCTION SUITE 

Remarque : Le mot-cle RESULT peut etre utilise dans toutes les fonctions pour definir le nom du 
resultat, qu'elles soient recursives ou non. Le type du resultat est toujours celui de la fonction. 

II.A.3 Les modules 

Un module est une unite de compilation qui a pour fonction de mettre a disposition, aux autres 
unites de compilation, des ressources logicielles comme des constantes, des variables, des types et des 
procedures. C'est une abstraction des donnees et des procedures qu'il encapsule, realisant ainsi une 
boite noire pour l'utilisateur et favorisant la maintenance des applications. Le module n'est pas 
executable, il n'y a que le programme principal qui le soit. La forme generale d'un module est la 
suivante : 

MODULE nom_du_module 

[instructions de specification] 
[CONTAINS 

[procedures globales]] 

[procedures internes]] 
END MODULE [ nom_du_module ] 
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Exemple II.5 

MODULE DONNEES 
IMPLICIT NONE 

REAL, PARAMETER :: Rgaz = 8.314e7 

CONTAINS 

SUBROUTINE INIT (rl, ul, Tl) 
IMPLICIT NONE 

REAL, INTENT (IN) :: rl, ul, Tl 
REAL : : gam, Mach 

gam = 5./3. 

! Mach = ul/Csound ; nombre de Mach 
Mach = ul * sqrt( 131/ (gam*Rgaz*Tl ) ) 
WRITE (6,*) "Le nombre de Mach est de ", Mach 

END SUBROUTINE INIT 
END MODULE DONNEES 

a) Acces aux ressources d'un module 

L'acces aux ressources d'un module se fait grace a l'instruction use nom_module. Cette ligne 
d'instruction se place toujours en tete de l'unite utilisatrice, avant le bloc declaratif et permet 
Fimportation de l'integralite du module dans cette unite de compilation. 

PROGRAM (ou SUBROUTINE, FUNCTION) 
USE Nom_de_Modulel 
USE Nom_de_Module2 
IMPLICIT NONE 

END PROGRAM (ou SUBROUTINE, FUNCTION) 

On peut toutefois restreindre l'acces a certaines ressources du module uniquement, en important 
partiellement le module. Pour ce faire, on place l'instruction only : suivie d'une liste exhaustive des 
ressources souhaitees, derriere la commande d'importation du module. 

USE Nom_de_Module, ONLY: varl, var2 , sous-programmel 

b) Renommer les importations 

Pour des raisons de lisibilite d'un programme, ou bien de compatibilite de plusieurs procedures, il 
est parfois souhaitable de renommer les entites importees. On utilise la double fleche qui est aussi le 
symbole d'affectation des pointeurs. 

USE Noml_de_Module => Nom2 , ONLY: varl => xy, var2 
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II.B FORMAT ET SYNTAXE 



II.B.l Le format 

En Fortran 90, on utilisera de preference le format libre, qui est une nouveaute de cette norme. Les 
lignes ont une longueur maximale de 132 caracteres. On peut coder plusieurs instructions sur une meme 
ligne en les separant par un point-virgule « ; ». Une instruction peut s'ecrire sur plusieurs lignes en 
utilisant le caractere de continuation esperluette « & ». Une chaine de caracteres ecrite sur plusieurs 
lignes doit etre precedee du « & » en plus de celui de continuation de la ligne precedente. 

Exemple II.6 

rr = rho_init ; uu = Cshock ; TT = Temp 
PRINT* , "Dans ce calcul le flux final est ", 
& Frad, "et l'opacite est ", 
& KRoss 

Exemple II.7 

PRINT* , "Dans ce calcul le flux final et& 
& l'opacite sont ", Frad, 
& KRoss 

II.B.2 Les identificateurs et la syntaxe 

Un identificateur est compose d'une suite de caracteres, limitee a 31, choisis parmi les lettres (non 
accentuees), de chiffres (sauf en premiere position) et du caractere « espace souligne » : « _ ». 
L'identificateur sert a nommer les objets que le programmateur veut manipuler. Le compilateur ne 
distingue pas les majuscules des minuscules, mais on peut toutefois s'en servir par soucis de lisibilite. 

Exemple II.8 : identificateurs autorises 
clight 

Nom_de_f ami 1 le 

TypeOfGas 

M2 

Exemple II.9 : identificateurs non autorises 
Norn de famille 
§TypeOfGas 
2M 

En Fortran 90, il n'existe pas de mots reserves, on peut done utiliser des mots cles du langage 
comme identificateur. Dans l'ideal, mieux vaut eviter de le faire, pour limiter les risques de confusion. 

Un espace est necessaire entre identificateurs, mots cles et nombres. 
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II.B.3 Commentaires 

Un commentaire commence par le caractere point d' exclamation « ! » et se termine par la fin de 
ligne. Done, il est possible d'inserer des commentaires apres des instructions, sur une meme ligne, mais 
pas avant. 

Exemple 11.10 

! gam : rapport des capacites de pression & volume 

! constants : Cp / Cv. 
gam = 5 . / 3 . 

! Mach = ul/Csound ; nombre de Mach 

Mach = ul * sqrt ( AtW/ ( gam*Rgaz*Tl ) ) ! adimensionne 

PRINT* , "Le nombre de Mach est de ", Mach 

gM2 = gam*Mach*Mach ! adimensionne 

II.B.4 Chaines de caracteres 

Les chaines de caracteres sont encadrees par des cotes simples « ' » ou doubles « " ». Elles 
neutralisent l'effet de tous les autres signes ou mots cles. 

Exemple II. 1 1 

! Le point d ' exclamation est utilise comme tel 
! dans une chaine de caracteres. dans ce cas 
PRINT* , "Le nombre de Mach est grand ! " 
"L' elephant est + grand que la souris." 
'3 * plus souvent ' 

n.C Les types et les declarations 

Comme toute information stockee en memoire est codee en binaire, il faut une cle pour connaitre 
la representation initiale de tel ou tel octet. Ainsi une chaine de caracteres, un nombre ou une valeur 
logique peuvent avoir dans l'absolu la meme representation binaire, suivant le codage utilise. C'est 
pourquoi, nous devons annoncer quel type d'objets nous manipulons, pour differencier leur codage 
binaire. Le Fortran 90 prevoit cinq types de bases et des combinaisons. 

logical, character, integer, real, double precision 
complex, type 

II est possible de definir de nouveaux types de donnees a partir d'un groupe de donnees appartenant 
a des types intrinseques, a l'aide de la declaration type. Nous verrons ce point ulterieurement. 

De plus on distingue les constantes, les variables scalaires, les vecteurs et tableaux au moyen 
d'attributs. Ces attributs servent a specifier ce que le programmateur souhaite faire de ces donnees 
declarees et permettent un controle de la bonne utilisation de ces donnees. 
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La forme generale de la declaration d'une variable est la suivante : 

type [ [ , liste_attributs ] :: ] liste_objets 

La liste des attributs est la suivante : 

parameter, dimension, allocatable, pointer, target, 
save, intent, optional, public, private, external, 
intrinsic 

Nous expliquerons leur signification au fur et a mesure de leur introduction dans les differents 
chapitre du cours. 

II.C.l Declaration de variables scalaires 

CHARACTER (LEN=12) :: elmt 
INTEGER, PRIVATE : : mmax, nmax 
INTEGER : : AtN ; REAL : : AtW 
REAL :: gM2 , CvM2 , ww, gam = 5./3. 
DOUBLE PRECISION :: x=5.d-4 
COMPLEX :: zl=(3,5), z2 

On voit que Ton peut initialiser les variables lors de leur declaration. Toute variable initialisee 
devient permanente, elle recoit done implicitement l'attribut « save ». 

II.C.2 Declaration de constantes 

Ces objets recoivent l'attribut « parameter », et doivent obligatoirement etre initialises lors de leur 
declaration. II devient alors interdit de changer leur valeur dans une instruction. 

CHARACTER (LEN=12) :: elmt= ' nomdef ami lie ' 

REAL, PARAMETER :: clight = 2.99791el0, Rgaz = 8.314e7 

COMPLEX, PARAMETER :: zl=(3,5), z2=(6,9) 

DOUBLE PRECISION, PARAMETER :: x=5.d-4 

Toutefois, certains programmateurs auront recours au typage par defaut. On commence par une 
declaration du type : 

IMPLICIT REAL (a-h, o-z) 

Ainsi tous les identificateurs commencant par une lettre de cet ensemble alphabetique (de a a h ou 
de o a z) seront automatiquement considered comme des reels. Sauf . si une declaration contraire est 
faite. Cette methode de declaration est pratique au commencement, car moins fastidieuse. Mais elle est 
source d'innombrables erreurs ! Par exemple, elle ne detecte pas les fautes de frappe, pourtant si 
nombreuses. Le Fortran 90 est justement un langage developpe dans le but de prevenir au maximum des 
erreurs de programmation, done il est absurde de se passer de cet atout. 

C'est pourquoi, il est vivement conseille d'inserer l'instruction : 
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IMPLICIT NONE 

avant tout bloc de declaration. Ce qui oblige le programmateur a declarer explicitement chaque variable. 

La declaration a l'aide du mot cle TYPE sert a creer des types derives definis par le 
programmateur. Nous verrons ce point ulterieurement, dans le chapitre II.D. 

II.C.3 Sous-type : la precision des nombres 

a) Avec le mot cle KIND 

Les types de base en Fortran 90 sont en fait des noms generiques renfermant chacun un certain 
nombre de sous-types que Ton peut selectionner a l'aide du mot cle KIND, lors de la declaration des 
objets. Pour chaque type predefini correspond un sous-type par defaut, c'est la simple precision en 
general. 

Les differentes valeurs du parametre KIND sont dependantes du systeme utilise et elles 
correspondent au nombre d'octets desires pour coder l'objet declare. Pour les chaines de caracteres, 
KIND sert a indiquer combien d'octets sont necessaires pour coder chaque caractere ! C'est utile dans 
des alphabets complexes. 

Dans une declaration, on indique soit directement le nombre d'octets, soit on donne un exemple du 
sous-type desire. 

Exemple 11.12 

REAL (KIND=8) :: x reel code sur 8 octets : 

double precision pour la plupart des calculateurs , 
et simple precision sur Cray C90, par exemple. 

INTEGER (KIND=2) :: k entier code sur 2 octets. 

En donnant le sous-type directement en exemple, au lieu du nombre d'octets. 

Exemple 11.13 

REAL, (KIND = KIND(l.dO)) :: x 

reel double precision quelque soit la machine utilisee. 

b) Avec des suffixes ou prefixes 

Une autre maniere de preciser le sous-type desire lors de l'ecriture des variables est de leur ajouter 
un suffixe pour les constantes numeriques, ou bien un prefixe pour les chaine de caracteres, qui contient 
la valeur du sous-type choisi. On insere le caractere « _ » entre la variable et le suffixe ou prefixe. 
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Exemple 11.14 

123456_2 entier code sur 2 octets 

34.567005678001_8 reel code sur 8 octets 

2_'fleur' chalne de caracteres, chacun code 

sur 2 octets 

ou bien 



Exemple 11.15 

INTEGER, PARAMETER : : 

123456_deux 

34. 567 00567800 l_huit 

deux 'fleur' 



deux = 2 , huit = 8 

entier code sur 2 octets 

reel code sur 8 octets 

chalne de caracteres, chacun code 

sur 2 octets 



II.C.4 Les fonctions intrinseques liees a la precision 



a) selected int kind (m) 

Cette fonction retourne un parametre de sous-type qui permet de representer les entiers n, en 
fonction de l'etendue choisie a l'aide de Tender m, tels que : 

-10 m < n < 10 m 
Elle retourne -1 si aucun parametre de sous-type ne repond a la demande. 



Exemple 11.16 

SELECTED_INT_KIND (30) 
SELECTED_INT_KIND (10) 
SELECTED_INT_KIND (10) 



retourne -1 
retourne 8 
retourne 4 



b) selected real kind (p. r) 

Cette fonction recoit au moins un des deux arguments p et m qui sont respectivement la precision 
et l'etendue choisies. Elle renvoie un parametre de sous-type qui permet de representer les reels 
repondant a cette precision, tels que : 

10" r < |x| < 10 r 

La fonction renvoie : 

- -1 si la precision demandee n'est pas disponible ; 

- -2 si l'etendue n'est pas disponible ; 

- -3 si la precision et l'etendue ne sont pas disponibles. 
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Exemple 11.17 



SELECTED_ 


_REAL_ 


KIND (2, 30) 


retourne 


4 


SELECTED_ 


_REAL_ 


_KIND(20,30) 


retourne 


16 


SELECTED_ 


_REAL_ 


_KIND(10, 10) 


retourne 


8 


SELECTED_ 


_REAL_ 


_KIND(40,10) 


retourne 


4 


SELECTED_ 


_REAL_ 


_KIND(20,5000) 


retourne 


-2 



c) range et precision 

En complement, l'etendue d'un nombre est connue a l'aide de la fonction RANGE. 

La precision d'un reel (nombre significatif de chiffres derriere la virgule) est connue avec la 
fonction PRECISION. 

Exemple 11.18 

RANGE (1234567890) retourne 9 

RANGE (123456. 7890) retourne 37 

PRECISION(l. 234567890) retourne 6 



II.D AUTRES TYPES 

Dans un souci de fiabilite dans l'execution et l'utilisation d'un programme, ou une partie, on est 
amene a definir un ensemble de variables. Le fait de les ranger ainsi sous un meme nom de type et dans 
un seul type previent d'erreurs de manipulation de ces variables et facilite leur declaration. Les tableaux 
sont deja des objets qui regroupent en ensemble de variables, mais c'est un cas particulier qui impose 
que toutes ces variables soient strictement du meme type. 

Dans le cas d'objets plus complexes, regroupant des donnees (champs ou composantes) de 
differents types, nous allons definir un type derive (ou structure) qui etend les types predefinis, et qui 
permet de manipuler ces objets globalement. Chaque composante est identified par son nom, elle peut 
etre soit de type intrinseque soit egalement de type derive. 

II.D.l Les tableaux 

Un tableau est un ensemble de variables qui sont toutes du meme type, il peut etre predefini ou de 
type derive. Chacune des variables est reperee par un indice numerique. La forme declarative generale 
est la suivante : 

TYPE, DIMENSION (...) [ , liste_attributs ] :: nom_tableau 

Un tableau peut avoir une taille statique, dans ce cas, sa dimension est explicitement ecrite lors de 
la declaration. Ou bien il a une taille dynamique, c'est-a-dire qu'elle sera definie en fonction de 
l'execution du programme, dans ce cas, la dimension est suggeree par le caractere deux-points « : ». 
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Dans le paragraphe consacre aux tableaux et a leur manipulation, nous reviendrons sur les formes 
declaratives des tableaux de taille dynamique. 

CHARACTER (LEN=8), DIMENSION (15) :: mois 
INTEGER, DIMENSION (12,31) :: jour 

LOGICAL, DIMENSION (2) :: reponse = (/. true .,. false ./ ) 

REAL, DIMENSION (:,:,:) :: position 

INTEGER, DIMENSION(4) , PRIVATE :: hmn 

REAL, DIMENSION (:), ALLOCATABLE :: TeV, rho 

II.D.2 Les types derives 

a) Construction d'un type derive (ou structure) 

Pour construire une nouvelle structure, on utilise la syntaxe suivante : 
type nom_type 

bloc_declaratif_des_champs 



end type nom_type 

Restriction du bloc declaratif : 

- L'attribut PARAMETER n'est pas autorise au niveau d'un champ. 

- L'attribut ALLOCATABLE est interdit au niveau d'un champ. 

- L'attribut POINTER est autorise au niveau d'un champ, mais pas TARGET. 

- L'initialisation d'un champ n'est possible qu'en Fortran 95. 

Attention : 

- Un objet de type derive est toujours considere comme un scalaire. meme si un champ possede 
l'attribut DIMENSION. 

- Les tableaux d'objets de type derives sont autorises et ils peuvent beneficier de l'attribut 
ALLOCATABLE. 

- Definition multiple de la structure : si le type derive n'est pas defini de maniere unique, par 
exemple lorsque la structure est passee en argument de procedure, l'attribut SEQUENCE est 
obligatoire. II assure un stockage en memoire dans le meme ordre des composantes en 
memoire. 

b) Manipulation d'un objet de type derive 

Pour utiliser des variables du type derive precedemment defini, on declare simplement : 
type (nom_type) [ , liste_attributs ] :: nom_var 

Pour acceder aux champs de l'objet que nous venons de definir, nous utiliserons le caractere pour 
cent « % ». Pour le reste des operations, la structure ainsi definie se manipule comme n'importe quelle 
variable. 
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c) Exemples 



Par exemple, un veterinaire voudrait creer un type derive « animal », dans lequel seraient recenses 
la race, Page, la couleur et l'etat de vaccination de ses « patients ». 

Exemple 11.19 

TYPE animal 

CHARACTER (LEN=2 0) :: race 

REAL : : age 

CHARACTER (LEN=15) :: couleur 
LOGICAL, DIMENSION(8) :: vacc 
END TYPE animal 

II peut ensuite manipuler globalement l'objet animal, ou bien acceder a chacun de ses champs, 
grace a l'operateur « % ». 

Pour declarer de nouveaux animaux, il utilise la forme suivante : 

Exemple 11.20 

TYPE (animal) :: ponpon, rox, roucky 

Pour rentrer les informations concernant ces nouveaux patients : 

Exemple 11.21 

ponpon%race = ' chat_de_gouttiere ' 

ponpon%age = 12 

ponpon%couleur = ' noir' 

ponpon%vacc = (/ .true. ,. false. , .true., .true.,& 

.false. ,. false., .true., .true./) 

II.D.3 Vers un langage oriente objet 

Fortran 90 n'est pas un langage oriente objet . a proprement parle, mais il s'en rapproche, car il 
possede des extensions objet qui accroissent la fiabilite des programmes et procurent un confort de 
programmation. En fait, il lui manque quelques notions, comme la hierarchie de types derives avec 
heritage. La norme Fortran 2000 les prevoit. 

Fortran 90 facilite cependant la manipulation des structures (objets de type derive), car il permet : 

- la definition d'autres fonctions de manipulation des structures, comme : 

- la surcharge des operateurs +,-,*, 

- de nouveaux operateurs en les associant aux fonctions correspondantes, 

- des autres fonctions ou procedures agissant sur ces structures. 

- d'encapsuler le type et les fonctions operant sur les objets de ce type, dans un module separe. 
Ce qui permettra d'y acceder facilement, avec l'instruction USE nom_module 

- de rendre privees les composantes interne d'une structure, avec l'instruction PRIVATE. 
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Ces differents point seront abordes ulterieurement. 

II.E LES fonctions intrinseques 

Fortran 90 possede une librairie tres complete de fonctions sur les scalaires. On trouvera bien sur 
toutes les fonctions mathematiques habituelles sur les entiers, les reels, et les complexes. Nous ne 
presenterons pas ces fonctions dans ce cours. 

De plus Fortran 90 possede des fonctions : 

- de manipulation de bits ; 

- de manipulation de chaines de caracteres ; 

- de conversion de type ; 

- d'interrogation sur les variables ; 

- des fonctions intrinseques liees a l'utilisation de tableaux (voir le Chapitre IV.C). 

Ces fonctions sont nommees et expliquees tout au long du cours, suivant le sujet developpe. 
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III Structures de controle 



A l'origine, les instructions d'un programme sont sequentielles : elles sont executees au fur et a 
mesure de leur ordre d' apparition. Mais rapidement, un programme plus elabore a besoin d'une 
souplesse de controle : dans tel cas faire ceci, par contre dans tels autre cas faire cela, le refaire n fois, 
etc. On a done besoin d' instructions pour faire des choix et des boucles. 

Le programme sera alors ecrit comme une suite de blocs. Un bloc peut-etre contenu dans un autre, 
mais seulement au niveau inferieur. II n'y a pas d'imbrication possible, sur le meme niveau de priorite. 
Ces suites de blocs forment le langage structure et augmentent la lisibilite du programme. 

Un bloc est compose d'une instruction initiale contenant le mot cle de la structure, de mots cles 
intermediaries et est ferme par un END mot cle de la structure. 

Les mots cles definissant les blocs sont : 

IF... THEN, ELSE IF... THEN, END IF 
DO, CYCLE, EXIT, END DO 
DO WHILE, END DO 

SELECT CASE, CASE, CASE DEFAULT, END SELECT 
WHERE, ELSEWHERE, END WHERE 
FORALL, END FORALL 



III.A LE CHOIX AVEC IF 

La construction typique d'un bloc est la suivante : 

IF (expression_logique) THEN 
blocl 

ELSE 

bloc2 
END IF 

Avec ses variantes : 

IF ( express ion_logiquel ) THEN 
blocl 

ELSE IF (expression_logique2) THEN 
bloc2 

END IF 
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Ou plus directement : 

IF (expression_logique) THEN 

blocl 
END IF 

Et encore plus simple, en une seule phrase : 

IF ( expression_logique ) blocl 

Le bloc IF peut etre etiquete, comme les autres structures presentees dans la suite. Dans ce cas, 
l'etiquette doit etre repetee au niveau du END, a sa suite. L'etiquette est un identifiant de la structure, 
c'est-a-dire un nom ; lorsqu'elle est optionnelle, elle augmente la lisibilite du programme ou bien elle 
est rendue obligatoire par la logique du programme. La syntaxe est la suivante : 

nom : IF (expression_logique) THEN 

blocl 
END IF nom 

Exemple III. 1 

stab = f(l,2) * f(2,2) 

IF (stab >= 1.01) THEN 

ecart = stab - 1. 

h = 0.5 * h 
ELSE IF (stab <= 0.99) THEN 

ecart = 1 . - stab 

h = 0.5 * h 

ELSE 

count = count + 1 
z = z + h 
f(3,l) = f(l,2) 
END IF 

III.B LE CHOIX AVEC SELECT CASE 

C'est une construction completement nouvelle en Fortran, elle permet de faire un choix multiple, 
sans expression logique, car le test porte sur une valeur unique de la variable. 

La syntaxe generale est la suivante : 
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[nom :] SELECT CASE (var) 
CASE (vail) 

blocl 
CASE (val2) 
bloc2 

CASE DEFAULT 

bloc default 
END SELECT [nom] 

Ou vail, val2.. sont des valeurs que peut prendre la variable selectionnee var. 

Le cas par defaut englobe tous ceux qui n'ont pas ete explicitement mentionnes, il est optionnel. 

Exemple III.2 

SELECT CASE (K) 

CASE (1) ! On integre en montant. 

PRINT* ," Integration directe : " 

h = Long* 1.2 / REAL ( nstep ) 

f(l,l) = 1. ; f(2,l) = 1.5 ; f(3,l) = 0. 

z = 0. 

nom_f ichierl = name 

CASE (2) ! On integre en descendant. 
PRINT* ," Integration inverse : " 
h = -Long* 1.2 / REAL (nstep) 
f(l,l) = r2/rl ; f(2,l) - u2/ul ; 

f(3,l) = T2/T1 

z = Long 

nom_f ichier2 = name 
END SELECT 
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Exemple III.3 

element : SELECT CASE (AtN) 

! AtN : n° atomique - AtW : Poids atomique 

CASE (1:9) 

elmt = ' trop_leger ' ; AtW = 0. 
CASE (12) 

elmt - 'magnesium' ; AtW = 24.305 
CASE (14) 

elmt = 'silicium' ; AtW = 34.855 
CASE (18) 

elmt = 'argon' ; AtW = 39.948 
CASE (47) 

elmt = 'argent' ; AtW = 107.868 
CASE (54) 

elmt = 'xenon' ; AtW = 131.2 9 
CASE DEFAULT 

elmt = ' non_utilise ' ; AtW = 0. 
END SELECT element 

IILC LE CHOIX AVEC WHERE 

Cette instruction ne s'applique qu'aux tableaux. L'instruction WHERE permet d'effectuer un 
choix, a l'aide d'un filtre, puis d'effectuer des instructions uniquement aux elements filtres. Ce bloc 
ressemble beaucoup a un bloc IF, mais il permet en plus une ecriture compacte pour les tableaux. 

La syntaxe generale est la suivante : 

WHERE (expression_logique) 

blocl 
ELSEWHERE 

bloc2 
END WHERE 

Remarque : Les blocs d'instructions blocl et blocl ne peuvent agir que sur des affectations 
concernant des tableaux. Dans le membre de droite des instructions, les tableaux sont obligatoirement 
conformant avec celui de l'expression logique. 

WHERE (tabl > 0.) 

tab2 = SQRT (tabl) 
ELSEWHERE 

tab2 - 0. 
END WHERE 
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En Fortran 95, la syntaxe d'un bloc WHERE s'etend grace a l'ajout de filtre de choix logique 
ELSEWHERE(filtre), de facon identique au bloc IF, ELSE IF, ELSE, END IF ; de plus il peut etre 
etiquete comme les structures precedentes. Le bloc WHERE s'ecrit alors sous la forme suivante : 

[nom : ] WHERE ( express ion_logiquel ) 
blocl 

ELSEWHERE (expression_logique2) 
bloc 2 

ELSEWHERE 
bloc3 
END WHERE [nom] 

III.D LABOUCLEDO 

Tres souvent en programmation, nous avons besoin de reproduire un grand nombre de fois les 
memes instructions. La boucle DO sert alors a iterer. 

La syntaxe generale est la suivante : 

[nom :] DO [bloc de controle] 
bloc 

END DO [ nom ] 

III.D. 1 Clause de controle iterative 

Si d'avance, on sait combien de fois, on a besoin de faire la meme serie d'operations, on utilise un 
compteur entier. 

[nom :] DO compt = intl, int2 [, int3 ] 
bloc 

END DO [nom] 

Ou intl est la valeur entiere de depart du compteur ; 

- intl est la valeur entiere d'arrivee ; 

- int3 est l'increment entier, par defaut il vaut 1 et est done optionnel ; 

- et compt est une variable entiere qui prendra les valeurs de intl a intl, en ajoutant int3 a 
chaque fois. Bien entendu, la derniere iteration ne tombe pas forcement sur la valeur de intl. 

Le nombre d'iterations est le plus grand nombre entre ( intl + int3 - intl ) I int3 et 0. 

Exemple III.4 

DO num = 1, 300 

T(num) = dx * num + y 

END DO 
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Attention : num = 4 en fin de boucle ! 

Exemple III.5 

alenvers : DO i = 201, 2, -1 
ux - (i - 1) * 5.e-3 

Tx = FUNTEMP (ul, ux, Tl, Tx, GT) 
f onction 

RT = FUNREST (ul, ux, Tl, Tx, GT) 
f onction 

PRINT* , ux, Tx, RT 

END DO alenvers 

III.D.2 Clause de controle WHILE 

Dans ce cas, on ne connait pas d'avance le nombre d'iterations, mais on connait plutot un domaine 
dans lequel on veut que la serie de calculs soit execute. Ce domaine est defini par une expression 
logique et 1' iteration continue tant que cette expression est VRAIE. 

Structure generale de la boucle DO WHILE : 

[nom :] DO WHILE (expression_logique) 
bloc 
END DO [nom] 

Exemple III.6 

Integre_choc : DO WHILE (n < nstep) 

rr = f(l,l) * rl ; TT = f(3,l) * Tl 

CALL RUNGEKUTTA (z, h, f, 3, n, Gas) 

n = z / h 
END DO Integre_choc 

Cette clause de controle est assez souvent utilisee pour lire les donnees d'un fichier dont on ne sait 
pas le nombre d'enregistrement. On utilise pour cela le marqueur de fin de fichier « eof », associe au 
parametre « iostat ». Nous verrons plus loin, avec les acces de fichiers et les entrees-sorties, comment 
on precede. Mais le boucle prend la forme suivante : 

Exemple III.7 

DO WHILE (eof == 0) 

READ (unit=65, iostat=eof), x, y, z 
END DO 



! appel de 
! appel de 
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III.D.3 Boucle infinie et alteration 

On peut aussi construire des boucles, qui sont a priori infinies, mais dans laquelle on insere une 
clause de sortie. Bien entendu, aucune boucle infinie n'est concevable, il faut l'interrompre par un 
moyen ou un autre. 

Pour cela on utilise 1' instruction EXIT : 

[nom :] DO [bloc de contrdle] 
bloc 

IF (expression_logique) EXIT [nom] 
END DO [ nom ] 

Pour alterer le deroulement d'une boucle, il est possible de forcer l'execution du programme a 
passer a la boucle suivante, avec l'instruction CYCLE. 

[nom :] DO [bloc de contrdle] 
blocl 

IF (expression_logique) CYCLE [nom] 
bloc2 

END DO [ nom ] 

III.E LE SCHEMA FORALL 

D'un point de vue algorithmique, le schema FORALL remplace certaines boucles DO. Mais pour 
le compilateur, il facilite la parallelisation de cette partie de code. La syntaxe du schema FORALL est la 
suivante : 

[nom :] FORALL (bloc de contrdle [, bloc de contrdle] 
[, f litre]) 

bloc 

END FORALL [nom] 

Le bloc de controle est une commande sur les indices d'un ou plusieurs tableaux. Si plusieurs 
plages d'indices sont definies a l'aide de plusieurs blocs de controle, le domaine doit etre rectangulaire. 
C'est-a-dire que chaque commande sur les indices ne doit pas dependre des autres indices. La 
construction des tableaux, ainsi que leur manipulation est detaillee dans le chapitre suivant. 

Exemple III.8 

FORALL (i=3:10, j=2:5) mat(i,j) - 12 

Exemple III.9 

FORALL (i=l:n, j=l:m, matl>0) 

mat2(i,j) = log (matl ( i , j ) ) 

END FORALL 
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IV Les tableaux 



Claire Michaut 



IV.A GESTION memoire des tableaux 

Un tableau est un ensemble d'elements ayant meme type et meme valeur de parametre de sous- 
type. II est defini par son type, son rang, ses etendues, son profil, sa taille. Un tableau peut contenir 
jusqu'a 7 dimensions. 

- Le rang {rank) d'un tableau est son nombre de dimensions. 

- L'etendue du tableau dans une dimension est le nombre d'elements dans cette dimension. 

- Le profil {shape) d'un tableau est un vecteur de rang 1 et d'etendue le rang du tableau et dont 
chaque element contient l'etendue dans la dimension correspondante. 

- La taille {size) du tableau est le produit des etendues sur toutes les dimensions. 

Deux tableaux de meme profil sont dits conformant. 

Dans le chapitre II.D, nous avons evoque la declaration des tableaux. Distinguons les tableaux qui 
ont une allocation statique de ceux qui ont une allocation dynamique. 

IV.A.l L' allocation statique 

II suffit de preciser le type du tableau, puis l'attribut dimension lors de la declaration en indiquant 
uniquement l'etendue de chaque dimension. 

type, DIMENSION ( b_inf 1 : b_supl , b_inf 2 : b_sup2 ...) :: nom_tab 

Ou b_infl, b_supl sont des entiers qui delimitent l'etendue de la premiere dimension du tableau 
nom_tab. II sont respectivement l'indice minimum et l'indice maximum. II en va de meme pour la 
seconde dimension avec les entiers b_inf2, b_sup2. 

IV.A.2 L' allocation dynamique 

On parle d'allocation dynamique lorsque Ton a pas defini au prealable l'etendue du tableau. II se 
peut, en effet, que cette donnee ne soit pas connue lors de l'ecriture du programme, par contre elle le 
deviendra pendant l'execution du programme. Ce phenomene est frequemment rencontre lorsqu'on lit 
des donnees dans un fichier genere par un autre programme. En Fortran 90, nous avons alors la 
possibilite de gerer de facon dynamique la memoire allouee a un tableau. II existe pour cela deux 
manieres distinctes : soit le passage en argument d'un tableau, et de sa dimension, dans une procedure, 
soit l'allocation directe de memoire. 
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La difference principale entre ces deux methodes est que : 

- dans la premiere, l'emplacement memoire des tableaux est alloue des l'entree dans la 
procedure et est libere a sa sortie ; 

- dans la seconde c'est le developpeur qui decide du moment de l'allocation et du moment de la 
liberation de la memoire. 

a) En argument de procedure ; les tableaux automatiques 

Ce cas se presente lorsque au moment du developpement d'une procedure, on ne connait pas 
encore la dimension du tableau, mais le programme appelant cette routine sera capable d'en fournir la 
valeur. Ou bien lorsque le meme sous-programme est appele a etre execute, avec des tableaux de tailles 
tres differentes. 

Avec passage en argument de la taille, la syntaxe generale se presente sous la forme suivante : 
SUBROUTINE nom (argl, nl, arg2 , n2 ...) 
IMPLICIT NONE 
INTEGER :: nl, n2 
type, DIMENSION (nl) :: argl 
type, DIMENSION (n2) :: arg2 

Sans passage de la taille en argument, la syntaxe generale se presente sous la forme suivante : 
SUBROUTINE nom (argl, arg2 , ...) 
IMPLICIT NONE 

type, DIMENSION ( : ) : : argl 
type, DIMENSION ( : ) : : arg2 

Pour connaitre l'etendue des tableaux, on a recourt a l'instruction SIZE. Mais, dans ce cas, 
l'utilisation d'une interface explicite est obligatoire, sinon le compilateur ne peut pas detecter des 
erreurs potentielles dans la manipulation de ces tableaux. Nous reviendrons sur ce point lorsque nous 
verrons la notion d'INTERFACE. 

b) Avec l'instruction ALLOCATE ; les tableaux dynamiques 

Pour generer dynamiquement des tableaux et liberer la memoire au moment opportun, Fortran 90 
possede l'instruction ALLOCATE et DEALLOCATE. Au prealable, le tableau considere devra etre 
declare avec l'attribut ALLOC ATABLE. 
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La forme generale a Failure suivante : 

type, DIMENSION( :[,...]) , ALLOCATABLE :: nom_tab 
INTEGER : : ok 

ALLOCATE ( nom_tab( diml [,...] ) [, stat = ok]) 

... [ bloc_instructions ] 

DEALLOCATE (nom_tab [, stat = ok]) 

L' utilisation du parametre stat est optionnel, c'est un controle du bon deroulement de 1' allocation 
ou de la liberation de la memoire. Lorsque l'allocation (ou la liberation) s'est correctement effectuee, 
l'instruction ALLOCATE (ou DEALLOCATE) renvoie pour la variable entiere ok la valeur 0, sinon ok 
prend une valeur superieure a 0. 

Exemple IV. 1 

SUBROUTINE SOLVER 

IMPLICIT NONE 

INTEGER : : npas 

REAL, ALLOCATABLE, DIMENSION ( : ) :: Prad, Frad 



ALLOCATE (Prad(npas), Frad (npas), stat = ok) 

IF (ok > 0) THEN 

PRINT*, "II y a un probleme d' allocation dans 
SOLVER" 

STOP 

END IF 



READ (18,*) Prad, Frad 

DEALLOCATE (Prad, Frad) 
END SUBROUTINE SOLVER 

II existe en plus une fonction qui permet de s'assurer que la memoire du tableau que Ton manipule 
est allouee ou non : 

allocated (array) 

Elle renvoie la valeur logique « .TRUE. » si le tableau allouable array est alloue, sinon la valeur 
« .FALSE. ». 
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IV.B LA MANIPULATION DES TABLEAUX 



a) Generalites 

Fortran 90 permet la manipulation globale des tableaux, ou de section de tableaux, en indiquant 
seulement leur nom. Cependant les tableaux doivent etre conformants pour etre utilises dans les memes 
instructions. Tous les operateurs utilises avec des variables scalaires sont alors utilisables avec des 
tableaux. 

L'initialisation d'un tableau au moment de sa declaration est permise a l'aide d'un c onstructeur de 
tableaux . C'est un vecteur de scalaires borne par les signes « (/ » et « /) ». 

Remarque importe : la valeur d'une expression faisant intervenir un tableau est entierement 
evaluee avant son affectation. 



b) Globale 

Attention toutefois a la simplicity de l'ecriture, car bien qu'elle allege les lignes de programme, il 
ne faut pas oublier que derriere un simple identifiant peut se cacher un tableau monstrueux ! D'un autre 
cote, cette ecriture compactee fait ressortir l'essentiel de la programmation et evite bien des erreurs sur 
les bornes et dimensions de tableaux. 

Exemple IV.2 

REAL, DIMENSIONS) :: tab 

tab - 100. 

Chacun des 4 elements du tableau tab prend la valeur 100. 
Exemple IV.3 

REAL, DIMENSION(4) :: tabl, tab2 , tab3 

tabl = 100. 
tab2 = 2 + tabl 

Chacun des 4 elements du tableau tab2 prend la valeur 102. 

Exemple IV.4 

INTEGER : : i 

REAL, DIMENSIONS) :: tabl, tab2= ( / 1 , 2 , 3 , 4 ) / ) 

tabl = (/10. ,15. ,20. ,25./) 
tab2 = tab2 * tabl 

Chacun des 4 elements du tableau tab2 prend respectivement la valeur 10 ; 30 ; 60 ; 100. 
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Exemple IV.5 

INTEGER : : i 

REAL, DIMENSIONS) :: tabl, tab2=(/ (i, i=l,4)/) 

tabl = 100. ; tab2 - tab2 * tabl 

Chacun des 4 elements du tableau tab2 prend respectivement la valeur 100 ; 200 ; 300 ; 400. 

On peut aussi repere par le caractere « : » l'etendue du tableau. La notation est alors plus lourde, 
mais permet au developpeur de ne pas oublier qu'il manipule un tableau a n dimension a ce moment. 
Pour reprendre l'exemple precedent, on peut ecrire sous la forme suivante. 

Exemple IV.6 

INTEGER : : i 

REAL, DIMENSION(4) :: tabl, tab2=(/ (i, i=l,4)/) 

tabl(:) = 100. ; tab2 ( : ) = tab2 ( : ) * tabl(:) 

Dans le cas de tableau a plusieurs dimension, on mentionnera chacune des dimensions par « : ». 
REAL, DIMENSION (4, 16) :: tabl, tab2 

tab2( : ) =20. / tabl( : , : ) 
c) Par section 

On peut aussi faire reference a une partie d'un tableau appelee section de tableau. Le sous-tableau 
ainsi defini est aussi un tableau. 

Les indices sont donnes soit par une suite arithmetique lorsqu'on a une section reguliere du tableau 
initial, soit par un vecteur d'indices, lorsqu'on veut extraire une section non reguliere. 

La section reguliere prend une forme typique : b_inf : b_sup : inc 

Ou b_inf et b_sup sont respectivement les indices inferieure et superieure de la section choisie, et 
inc est un entier qui represente l'increment. Par defaut inc vaut 1. C'est le meme principe que dans la 
boucle DO avec clause iterative. 

Exemple IV.7 

INTEGER : : i 

REAL, DIMENSION (20) :: tabl = (/ (i, i=l,20)/) 
REAL, DIMENSION(5) :: 

tab2( : ) = tabl (1:20:4) 
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De plus, on peut combiner les differentes manipulations possibles, en voici un certain nombre 
d'exemples. 

Exemple IV.8 

! on selectionne les elements de la ligne 2 
! et des colonnes de 4 a 8 

tabl(2, 4:8) 

! on selectionne tous les elements de la ligne 2 
tabl(2, : ) 

! on selectionne tous les elements des lignes 1, 2 et 3 
tabl( :3, : ) 

! on selectionne les elements des lignes 1, 3 et 5 
! et des colonnes de 4 a n 

tabl( 1:6:2,4: n) 

La section quelconque prend la forme d'un vecteur d'indices : (/il, i2, i3.../), a l'aide du 
constructeur de vecteurs. 

Exemple IV.9 

! on selectionne les elements des lignes 1, 5 et 19 
! et des colonnes de 4, 5 et 10 

tabl( (/1,5,19/), (/4,5,10/)) 

On peut aussi affecter directement une section avec l'identifiant d'un vecteur. 
Exemple IV. 10 

INTEGER, DIMENSIONS) :: vect_ind - (/3,6,7,9/) 
tabl(vect_ind, (/4,5,10/)) 

Lorsqu'un sous-tableau est constitue d'elements repetes, il ne peut pas recevoir une nouvelle 
affectation, done il ne peut etre place a gauche d'une expression. 



Exemple IV. 1 1 



INTEGER, DIMENSION(4) 
INTEGER, DIMENSION ( 10 ) 
INTEGER, DIMENSION(4) 
! on peut ecrire : 

tab2 = tabl ( vect_ind) 
! mais il est interdit d' ecrire : 
tabl ( vect_ind) = ... 



: vect_ind = (/3,6,3,9/) 
: tabl 
: tab2 
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IV.C LES fonctions intrinseques sur les tableaux 



Fortran 90 propose toute une gamme de fonctions intrinseques qui offre l'avantage de simplifier 
l'utilisation des tableaux. Nous distinguerons les fonctions passives qui renvoient seulement une 
information sur le tableau des fonctions actives qui modifient le tableau. 

Dans la suite, on utilisera comme argument les mots suivants : 

- array qui designe le tableau sur lequel on applique la fonction ; 

- vect qui designe un vecteur ; 

- mask qui designe un tableau d' expressions logiques conformant au tableau array ; 

- dim qui restreint la fonction a s'appliquer a la dimension mentionnee. 

IV.C.l Les fonctions d' interrogation 

maxloc (array [, mask] [, dim]) 
minloc (array [, mask] [, dim]) 

Ces fonctions retournent respectivement un tableau de rang 1 correspondant a la position de 
l'element maximum et de l'element minimum du tableau array. Fortran 95 admet en plus la precision 
de la dimension dim, ce que ne fait pas Fortran 90. 



Ces fonctions retournent respectivement un tableau d'entiers du type par defaut et de rang 1 qui 
contient la liste des bornes inferieures du tableau array (ou de la dimension dim mentionnee) et la liste 
des bornes superieures du tableau array (ou de la dimension dim mentionnee). 

size (array [, dim]) 

Cette fonction retourne la taille totale du tableau array ou l'etendue le long de la dimension dim. 

shape (array) 

Cette fonction retourne un tableau d'entiers de rang 1 qui contient le profil du tableau array. 
Exemple IV. 12 



lbound (array [, dim]) 
ubound (array [, dim]) 



INTEGER, DIMENSION (-1:1, 3) 



tabl 



/-15 2 35 



tab\ = 



4 



50 -36 



V -7 28 -9 ; 
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maxloc (tabl) donne (/2, 2/) 

donne (/2, 3/) 

donne (/-l, 1/) 

donne (/l, 3/) 

donne 9 

donne (/3, 3/) 



minloc (tabl) 
lbound (tabl) 
ubound (tabl) 
size (tabl) 
shape (tabl ) 



allocated (array) 

Cette fonction fournit la valeur « .TRUE. » si le tableau array, declare avec allocatable, est alloue 
et la valeur « .FALSE. » sinon. 



IV.C.2 Les fonctions de reduction 

Les fonctions de reduction renvoient soit un scalaire soit un tableau d'un rang inferieur a celui 
passe en argument. 



a) Les fonctions all et any 
all (mask [ , dim] ) 

Cette fonction renvoie la valeur logique « .TRUE. », si tous les elements, pris un par un, repondent 
aux criteres du masque, sinon la valeur « .FALSE. ». Si dim est precise, la fonction renvoie un tableau 
qui contient les valeurs logiques « .TRUE. » ou « .FALSE. » suivant que les elements de cette 
dimension repondent aux criteres du masque. 

Exemple IV. 13 

tabl(l,:) = (/l, 3, 5, 7/) ; tab2 ( 1 , : ) = (/l, 2, 3, 4/) 
tabl(2,:) = {12, 4, 6, 8/) ; tab2 ( 2 , : ) = (/5, 6, 7, 8/) 
! reduction globale 

ALL (tabl==tab2) donne .false. 

ALL (tabl/=tab2) donne .false. 

! reduction par colonne 
ALL (tabl/=tab2, dim = 1) 

donne /(.false., .true., .true., .false./) 
! reduction par ligne 
ALL (tabl/=tab2, dim = 2) 

donne /(.false., .false./) 

Si on a deja des tableaux de variables logiques, on les utiliser directement dans l'instruction ALL. 
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Exemple IV. 14 

tabl(l,:) = (/.true., .true., .true., .true./) 
ALL (tabl) donne .true. 

On peut aussi tester les valeurs d'un tableau. 

Exemple IV. 15 

REAL, DIMENSION (5,10) :: tabl, tab2 

IF (ALL (tabl) >= 0) THEN 

tab2 = SQRT (tabl) 
END IF 

De facon symetrique a la fonction all, on peut utiliser la fonction suivante : 
any (mask [ , dim] ) 

Cette fonction renvoie la valeur logique « .TRUE. », si au moins un des elements, pris un par un, 
repondent aux cri teres du masque, sinon la valeur « .FALSE. ». Si dim est precise, la fonction renvoie 
un tableau qui contient les valeurs logiques « .TRUE. » ou « .FALSE. » suivant que les elements de 
cette dimension repondent aux criteres du masque. 

Exemple IV. 16 

tabl(l,:) = (/l, 3, 5, 7/) ; tab2 ( 1 , : ) = (/l, 2, 3, 7/) 
tabl(2,:) = {12, 4, 6, 8/) ; tab2 ( 2 , : ) = (/5, 6, 4, 8/) 
! reduction globale 

ANY ( tabl==tab2 ) donne .true. 

ALL (tabl/=tab2) donne .true. 

! reduction par colonne 
ALL (tabl/=tab2, dim = 1) 

donne /(.true., .true., .true., .false./) 
! reduction par ligne 
ALL (tabl/=tab2, dim = 2) 

donne /(.true., .true./) 

Si on a deja des tableaux de variables logiques, on les utiliser directement dans l'instruction ALL. 
Exemple IV. 17 

tabl(l,:) = (/.true., .true., .false., .true./) 
ANY (tabl) donne .true. 

On peut aussi tester les valeurs d'un tableau. 
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Exemple IV. 18 

REAL, DIMENSION (5,10) :: tabl, tab2 

IF (ANY (tabl) <= 0) THEN 

PRINT*, 'On ne peut pas calculer la racine carree. ' 

EXIT 
END IF 

b) La fonction count 
count (mask [, dim]) 

Cette fonction compte le nombre d'occurrences qui prennent la valeur logique « .TRUE. », suivant 
les criteres du masque. Si dim est precise, la fonction renvoie un tableau d'entiers contenant le nombre 
d'occurrences VRAI. 

Exemple IV. 19 

tabl(l,:) = (/l, 3, 5, 7/) ; tab2 ( 1 , : ) = (/l, 2, 3, 7/) 
tabl(2,:) = {12, 4, 6, 8/) ; tab2 ( 2 , : ) = (/5, 6, 4, 8/) 
! compte global 

COUNT (tabl==tab2) donne 3 

COUNT (tabl/=tab2) donne 5 
! compte par colonne 

COUNT (tabl/=tab2, dim = 1) donne /(l, 2, 2, 0/) 
! compte par ligne 

COUNT (tabl/=tab2, dim = 2) donne /(2, 3/) 

Si on a deja des tableaux de variables logiques, on les utiliser directement dans l'instruction ALL. 
Exemple IV.20 

tabl(l,:) = (/.true., .true., .false., .true./) 
COUNT (tabl) donne 3 

c) Les fonctions maxval et minval 
maxval (array [, dim] [, mask]) 
minval (array [, dim] [, mask]) 

Ces fonctions renvoient respectivement la valeur maximum et la valeur minimum du tableau 
array, ou le tableau des maxima et des minima suivant la dimension du tableau imposee et selon les 
criteres du masque s'il est precise. 

Note : Si le tableau est de taille nulle, le resultat est respectivement le plus grand nombre positif et 
negatif disponible sur l'ordinateur. 
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Exemple IV.21 

tabl(l,:) = (/l, 3, 5, 7/) ; tab2 ( 1 , : ) = (/l, 2, 3, 7/) 
tabl(2,:) = {12, 4, 6, 8/) ; tab2 ( 2 , : ) = (/5, 6, 4, 8/) 

! reduction globale 

MAXVAL (tabl) donne 8 

! reduction partielle, par colonne 

MINVAL (tabl,DIM=l,MASK=tabl>l) donne (/2, 3, 5, 7/) 

d) Les fonctions sum et product 
sum (array [, dim] [, mask]) 
product (array [, dim] [, mask]) 

Ces fonctions calculent respectivement la somme et le produit des elements d'un tableau de 
nombres ; ou suivant la dimension indiquee et selon les criteres du masque, s'il est precise. 

Note : Si le tableau est de taille nulle, le resultat est respectivement et 1 . 
Exemple IV.22 





'12 3 






HI) et tabl= 




3 






v 5 6 4 




SUM(tabl) 


donne 


36 




PRODUCT (tabl, MASK=tabl<6) 


donne 


120 




SUM(tab2, DIM=2) 


donne 


(/13 


, 23/) 


PRODUCT ( tab2 , DIM=2, MASK=tab2>4) 


donne 


</7, 


240/) 



IV.C.3 Les fonctions de construction et de transformation 



a) La fonction merge 
merge ( array 1, array2, mask) 

Cette fonction fabrique un tableau, a partir des deux tableaux conformants et de meme type arrayl 
et arrayl, en fonction des criteres imposes par le masque. Si la valeur du masque est vrai, on utilise 
l'element de arrayl, si la valeur du masque est faux on utilise l'element de array2. 

Exemple IV.23 





(1 2 


3^ 




1 10 20 


30 N 


tab\ = 


4 5 


6 


et tabl = 


40 50 


60 




J 8 


9j 




V70 80 


90 y 



MERGE (tabl, tab2, tabl>4) 



/10 20 30\ 



donne tab\ = 



40 5 6 
7 8 9 
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b) La fonction reshape 
reshape (array, shape [, pad] [, order])) 

Cette fonction construit un tableau, a partir des elements du tableau array, et dont le profil est 
donne par le tableau d'entiers positifs de rang 1 shape. L'attribut optionnel pad est un tableau du meme 
type de array, dont les elements seront utilises si la taille de shape est superieure a celle de array. 

Exemple IV.24 

RESHAPE ( (/(i, i=l,12)/), (/3, 4/)) 

(l 4 7 10 N 



donne 



2 5 8 11 

3 6 9 12 



c) Les fonctions pack et unpack 
pack (array, mask [, vector]) 

Cette fonction construit un tableau de rang 1, en compressant les valeurs du tableau array, suivant 
les criteres du masque. Si le vecteur optionnel vector est present, le resultat aura la taille de vector. 

unpack (vector, mask [, field]) 

Cette fonction est assez symetrique de la precedente. Elle construit un tableau de rang superieur a 
1 , en decompressant les valeurs du vecteur vector, suivant les criteres du masque. Si le tableau optionnel 
field est present, les valeurs de field, qui sont du meme type que celle de vector, seront utilisees lorsque 
le masque contient une valeur fausse. Le resultat aura le profil du masque et le type du vecteur. 

Exemple IV.25 

(l 4 7 10 x 



On a tab\= 



2 5 8 11 

3 6 9 12 



vectl = (/ (i, i=-12,-l) /) 

PACK (tabl, MASK=tabl<10, VECTOR=vectl ) 

donne (/l, 2, 3, 4, 5, 6, 7, 8, 9, -3, -2, -1 /) 

d) La fonction spread 
spread (array, dim, n) 

Cette fonction cree un nouveau tableau en dupliquant n fois les valeurs du tableau array, le long de 
la dimension mentionnee. Le resultat est un tableau d'un rang superieur a array. 

e) Les fonctions cshift et eoshift 
cshift (array, shift [, dim]) 

Cette fonction opere un decalage circulaire sur les elements du tableau array, le long de la 
dimension indiquee, d'autant d'indices que shift indique. 
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Attention : shift est un entier positif, pour un decalage vers les indices decroissants (on effectue i = 
i-shift), ou un entier negatifs, pour un decalage vers les indices croissants. Par defaut dim est a 1. 



eoshift (array, shift [, boundary] [, dim]) 

Cette fonction agit comme la precedente, mais donne en plus la possibilite de remplacer les 
elements perdus a la frontiere par des elements de remplissage donnes par le tableau boundary. 



Note : Si boundary n'est pas present, les valeurs de remplissage par defaut dependent du type du 
tableau array : 



- pour un INTEGER, la valeur par defaut est 

- pour un REAL, 

- pour un COMPLEX, 

- pour un LOGICAL, 

- pour un CHARACTER, 



0; 
0.0 ; 

(0.0,0.0) ; 
.FALSE. ; 

un chaine composee de blancs. 



Exemple IV.26 



CSHIFT (tabl, SHIFT=2, DIM=1) 



EOSHIFT (tabl, SHIFT=-1, DIM=2 ) 
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donne 
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3 
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9 



IV.C.4 Les fonctions propres aux matrices et vecteurs 

a) Multiplication de vecteurs et de matrices 

Pour effectuer le produit scalaire de deux vecteurs : 
dot_product (vectl, vect2 ) 

et la multiplication de deux matrices : 
matmul (mata, matb) 
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b) matrice transposee 

transpose (matrix) 

Cette fonction renvoie la matrice transposee du tableau matrix. 
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V 



Les pointeurs 



V.A 



INTRODUCTION 



V.A.1 



Definition : 



Le pointeur est un objet qui peut designer des objets differents au cours de l'execution d'un 
programme, il recoit l'attribut POINTER. II agit comme un alias et est du meme type que les objets 
qu'il pointe. L'objet designe par le pointeur est appele la cible du pointeur, il doit recevoir l'attribut 
TARGET. Pour assigner un objet a un pointeur, on utilise le symbole de la double fleche materialisee 
pas les caracteres : « => ». Un pointeur a aussi le droit d'etre valoriser vers un autre pointeur. 

La notion de pointeur est surtout utile dans le contexte de 1' allocation dynamique de memoire. 

V.A.2 Declaration d'un pointeur et d'une cible 

REAL, TARGET : : targ 
REAL, POINTER :: ptr 

Le type, le parametre de type et le rang du pointeur et de la cible doivent etre identiques. 

Et pour affecter targ comme cible de ptr (ou pour valoriser le pointeur) : ptr => targ. 

Remarque : pour declarer un tableau de pointeurs, seul son rang (le nombre de dimensions) doit 
etre declare et ses bornes prendront les valeurs de l'objet pointe. 



Exemple V. 1 



REAL, POINTER :: ptrl 

INTEGER, DIMENSION^,:), POINTER :: ptr2 

! attention ptr3 n'est pas un tableau de pointeurs. 

COMPLEX, POINTER :: ptr3 ( : ) 



V.B 



ETAT D'UN POINTEUR 



Un pointeur se trouve forcement dans un des trois etats suivants : 

- indefini : comme lors de sa declaration, 

- nul : il ne pointe sur rien, c'est-a-dire qu'il n'est pas l'alias d'une variable, 

- ou associe : il pointe sur une variable. 
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Note : en Fortran 95, on peut forcer un pointeur a l'etat nul, lors de sa declaration. 
REAL, POINTER :: ptr => NULL ( ) 



Claire Michaut 



a) La fonction nullify 

L'instruction NULLIFY permet de forcer un pointeur a l'etat nul. 
syntaxe : NULLIFY (ptrl) 

! norme 95 : remplacee par la valorisation 
ptrl => NULL ( ) 

b) La fonction associated 

L'etat d'un pointeur peut etre connu grace a la fonction intrinseque ASSOCIATED, dont la 
syntaxe generale est la suivante : 

associated (pointer [, pointer] [, target]) 

associated (pointer) 

Cette fonction renvoie « .TRUE. », si pointer est associe a une cible, s'il est dans l'etat nul, la 
fonction renvoie « FALSE. ». 

associated (pointerl, pointer2) 

Cette fonction renvoie « .TRUE. », si pointerl et pointerl sont associes a la meme cible, sinon elle 
renvoie « FALSE. ». Remarque : si pointerl et pointerl sont tous les deux dans l'etat nul, elle renvoie 
« FALSE. ». 

associated (pointer, target) 

Cette fonction renvoie « .TRUE. », si pointer est associe a la cible target, sinon elle renvoie 
« FALSE. ». 

Attention : le pointeur ne doit pas etre dans l'etat indetermine pour utiliser cette fonction ! 

c) Operation sur les pointeurs 

On a vu que l'operateur « => » sert a valoriser le pointeur. Le pointeur prend alors la valeur de la 

cible. 

L'operateur « = » sert a affecter une valeur, mais attention lorsque les operandes sont des 
pointeurs, 1' affectation porte sur la cible et non sur le pointeur. 
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L' instruction ALLOCATE, que nous avons deja vu pour 1' allocation dynamique de memoire des 
tableaux, sert aussi a associer un pointeur et a lui allouer de la memoire. L'instruction DEALLOCATE 
sert a liberer la memoire, comme pour un tableau. 

Notes : 

- Dans ce cas, l'espace alloue n'a pas de nom et la seule facon d'y acceder est d'utiliser le 
pointeur. Le pointeur est 1' alias direct du chemin qui mene a cette memoire. 

- Apres l'instruction DEALLOCATE, le pointeur est dans l'etat nul . 

- On ne peut pas executer l'instruction DEALLOCATE sur un pointeur dont l'etat est 
indetermine. 

- On peut aussi faire de 1' allocation dynamique de memoire pour des variables, par 
1' intermediate de pointeurs. 

Exemple V.2 : tri de chaine de caracteres 
MODULE CHAINE 

END MODULE CHAINE 
PROGRAM TRI_CHAINE 
IMPLICIT NONE 
CHARACTER (LEN=80 ) , 
CHARACTER (LEN=80 ) , 
CHARACTER (LEN=80 ) , 
INTEGER : : num, i 
LOGICAL : : f ini 

PRINT*, "Quel est le nombre de chaines ? " 
READ ( * , * ) num 

ALLOCATE ( chaine ( num ) , pt r_ch ( num ) ) 
lis : DO i = 1, num 

PRINT*, "Donnez votre message : " 

READ ( * , * ) chaine(i) 

ptr_ch(i) => chaine(i) 
END DO 



DIMENSION( : ) , TARGET :: chaine 
DIMENSION ( : ) , POINTER :: ptr_ch 
POINTER : : ptr_int 
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tri_tous : DO 

fini = .TRUE. 

tri_chacun : DO i = 1, num-1 

IF (chaine(i) > chaine(i+l)) THEN 
fini = .FALSE. 
ptr_int => ptr_ch(i) 

ptr_ch(i) => ptr_ch(i+l) 
ptr_ch(i+l) => ptr_int 
END IF 
END DO tri_chacun 
IF (fini) EXIT 
END DO tri_tous 
END PROGRAM TRI_CHAINE 

Dans l'exemple precedent, on aurait pu faire le meme tri a l'aide d'un tableau ALLOCATABLE, a 
la place du tableau de pointeurs. Neanmoins, il parait que le temps d'execution d'un programme est 
ameliore dans le cas de tableaux de grandes tailles et de longues chaines de caracteres ! Car la gestion 
de la memoire differe dans ces deux cas. 

V.C.I En argument de procedure 

II est possible de passer en argument de procedures un pointeur. Nous distinguons le cas ou le 
pointeur est deja dans la definition de la procedure, du cas ou l'argument muet n'est pas un pointeur. 

a) L'argument muet a l'attribut pointer 

Dans ce cas, le pointeur figure obligatoirement dans l'appel de la procedure et n'est pas forcement 
associe au moment de l'appel. L'information reelle qui passe dans la procedure est le descripteur du 
pointeur : c'est-a-dire son adresse, ses dimensions... 

Remarques : 

- L'interface doit etre explicite, pour que le compilateur sache que la procedure possede en 
argument muet un pointeur. 

- L'attribut intent ne peut pas etre utilise pour qualifier le pointeur qui est en argument muet. 

b) L'argument muet n'a pas l'attribut pointer 

Dans ce cas, le pointeur est obligatoirement associe avant l'appel de la procedure. C'est alors 
l'information concernant la cible qui passe dans la procedure, c'est-a-dire son adresse. On se retrouve 
alors dans un cas habituel de passage de variables dans un appel de procedure. 
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c) Cible en argument de procedure 

Si un argument en appel de procedure a l'attribut TARGET, tout pointeur l'ayant pour cible au 
moment de 1' appel, devient indefini au retour de la procedure. De facon generale, la norme ne garantit 
pas la conservation de l'association du pointeur entre sa cible passee en argument d'appel de procedure 
et la cible correspondante en argument muet. 
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VI CONTROLE DE VISIBILITE 



VI.A LES ressources privees et publiques 

II est souvent souhaitable que l'integralite des ressources d'un module ne soit pas accessible par 
les autres unites de programme. En effet, lorsqu'on definit un module, on souhaite faire partager un 
certain nombre de variables, de definitions de type, de procedures, via l'instruction use nom_module ; 
mais pour des besoins specifiques a ce module, on est en droit de definir des ressources propres a ce 
module. 

- Les ressources non visibles depuis l'exterieur sont dites privees (private). 

- Tandis que les ressources visibles depuis l'exterieur sont dites publiques (public). 

Par defaut. toutes les ressources d'un module sont publiques . 

On a souvent interet a rendre privees certaines ressources d'un module, cela evite les risques de 
conflits avec les ressources d'un autre module. 

Une autre facon de restreindre l'utilisation des ressources d'un module est d'inserer l'instruction 
only nom_var, lors de l'acces a un module, dans l'instruction use nomjnodule (voir plus haut 
chapitre II.A.3a)). 

a) Les instructions PUBLIC et PRIVATE 

Le mode par defaut est le mode public. Les instructions PUBLIC et PRIVATE, sans argument, 
permettent de changer le mode. Un telle instruction ne peut apparaitre qu'une seule fois dans un 
module et affecte l'ensemble du module. 

MODULE PARAM 
IMPLICIT NONE 

INTEGER, DIMENSION ( : ) :: par 

PRIVATE ! tout le module devient prive. 

REAL : : x, y 

INTEGER : : i, j 

END MODULE PARAM 

b) Les attributs PUBLIC et PRIVATE 

Si Ton veut rendre public ou prive n'importe quel objet d'un module (ou d'une procedure), on lui 
affecte soit l'attribut PUBLIC ou soit l'attribut PRIVATE, directement lors de sa declaration. 
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Syntaxe generate : 

type, private :: var 
type, public :: var 

Attention toutefois a ne pas confondre le fonctionnement des attributs PRIVATE et PUBLIC, avec 
celui des instructions PRIVATE et PUBLIC, malgre leur similitude. 

MODULE INIT 
IMPLICIT NONE 

INTEGER, DIMENSION ( : ) :: par 

PRIVATE ! rend 1 ' environnement prive . 

REAL : : x, y 
INTEGER : : i, j 

REAL, DIMENSION( : ) , PUBLIC :: vect 

! rend public le tableau "vect", done exportable. 
PUBLIC : : lect 

CONTAINS 

SUBROUTINE lect (x,y) 

! cette procedure est publique. 

END SUBROUTINE lect 
FUNCTION diet (vect) 

! cette procedure est privee. 

END FUNCTION diet 
END MODULE INIT 

VI.A.2 les types derives semi-prives 

On parle de type derive semi-prive lorsqu'un type derive, defini au chapitre II.D.2, est public, mais 
dont toutes les composantes sont privees. Son interet est de permettre au developpeur de modifier la 
structure du type derive, sans affecter les autres unites de programme qui l'utilisent. 

Les attributs PUBLIC et PRIVATE s'appliquent aux types derives, comme pour les autres types de 
variables. 

Alors un type derive est soit : 

- transparent : il est public et ses composantes sont aussi publiques ; 

- semi-prive : il est public et ses composantes sont toutes privees ; 

- ou prive. 
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Exemple VI. 1 

MODULE BASE 

TYPE STAR 

PRIVATE 

INTEGER : : num 

REAL, DIMENSION( : ) , POINTER :: lambda 
END TYPE STAR 
END MODULE BASE 

Dans cet exemple, de definition d'un type semi-prive, seules des procedures internes au module 
base peuvent acceder aux champs de la variable de type star. 

Exemple VI.2 

MODULE BASE 

TYPE STAR 

PRIVATE 

INTEGER : : num 

REAL, DIMENSION( : ) , POINTER :: lambda 
END TYPE STAR 
CONTAINS 

FUNCTION VALOR (bl, b2 ) 

IMPLICIT NONE 

TYPE (STAR) :: bl, b2 

bl%num = b2%num - 1 

valor = (bl% lambda + b2% lambda) / 2. 
END FUNCTION VALOR 
END MODULE BASE 

VI.A.3 Restriction de 1' instruction USE 

De la meme maniere que Ton peut rendre prive certaines ressources d'un module, on peut limiter 
Faeces aux ressources d'un module, lors de l'instruction USE, dans l'unite utilisatrice. 

Syntaxe generate : 

USE nom_module , ONLY : list_var 
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Exemple VI.3 

MODULE PARAM 

IMPLICIT NONE 

INTEGER, DIMENSION ( : ) :: par 
REAL : : x, y 
INTEGER : : i, j 
END MODULE PARAM 

PROGRAM MACHINE 

USE PARAM, ONLY : x, y 

! seules ces 2 variables sont accessibles 

END PROGRAM MACHINE 

II est aussi possible de renommer des ressources d'un module, si les noms proposes par le module 
ne conviennent pas a 1' unite utilisatrice : par exemple lorsqu'un de ces noms est deja affecte ! Pour 
renommer, on utilise l'operateur d'affectation « => », identique a celui des pointeurs. 

Syntaxe generate : 

USE nom_module, ONLY : old_var=>new_var, old_vec->new_vec 

Exemple VI.4 

MODULE PARAM 
IMPLICIT NONE 

INTEGER, DIMENSION ( : ) :: par 
REAL : : x, y 
INTEGER : : i, j 
END MODULE PARAM 

PROGRAM MACHINE 

USE PARAM, ONLY : x=>varx, y=>vary 
END PROGRAM MACHINE 
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VII CONTROLE DE COHERENCE 



VILA LES arguments de procedure 
VII.A.l Vocation des arguments 

Lors de la declaration des arguments de procedure, il est possible d'affiner le controle sur la 
manipulation de ces arguments, au cours de la procedure, en leur specifiant une vocation : 

- etre un parametre d'entree seulement, done ne pas recevoir de nouvelle affectation ; 

- etre un parametre de sortie seulement, done recevoir certainement une affectation ; 

- ou bien, etre a la fois un parametre d'entree et de sortie de la procedure. 

Pour specifier cette vocation a etre un parametre d'entree ou de sortie de procedure, on utilise 
l'attribut INTENT, qui peut etre in, out ou inout. 

Syntaxe generale : 

procedure nom_proc (argl, arg2 , argn) 

type, INTENT (in) :: list_argl 
(et/ou) 

type, INTENT (out) :: list_arg2 
(et/ou) 

type, INTENT (inout) :: list_arg3 

Ou list_argl, Ust_arg2 et Ust_arg3 sont des listes de variables prises parmi les arguments de la 
procedure. 



Exemple VII. 1 

SUBROUTINE INIT (rl, ul, Tl) 

IMPLICIT NONE 

REAL, INTENT (IN) :: rl, ul, Tl 
REAL : : Mach, Miso, Mrad, Mcon 
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! Calcul des constantes du systeme d' equations. 

! gam = rapport des gaz par f aits gam = Cp / Cv. 
gam = 5 . / 3 . 

! Mach = Cshock / Csound = ul/Cs ; nombre de Mach 

! adimensionne 

Mach = ul * sqrt (AtW/ ( gam*Rgaz*Tdeg*Tl ) ) 
WRITE (6,*) "Le nombre de Mach est de ", Mach 
gM2 = gam*Mach*Mach ! adimensionne 

CvM2 = 1.5 / (Mach*Mach) ! adimensionne 

! On a pose dans le calcul : w = (a Tl**4) / (3 rhol ul**2) 

ww = 7.568e-15 * (Tl*Tdeg)**4 / ( 3 . *rl*ul*ul ) ! 
[ erg/cm3 ] 

! On a aussi pose : 1 = c/ul / ( Krossland*L ) , mais on en a 
pas encore besoin. 

Miso = sqrt ( ( 3 . *gam-l . ) / ( gam* ( 3 . -gam) ) ) 

PRINT*, "Le seuil isotherme est ", Miso 

Mrad = 7**(7./6.) * ( 6*gam) ** ( -0 . 5 ) * (ww * gM2)**(- 
1./6. ) 

PRINT*, "Le seuil radiatif est ", Mrad 
Mcon = 2.38 * Mrad 

PRINT*, "Le seuil continu est ", Mcon 

END SUBROUTINE INIT 

Dans l'exemple precedent, les arguments rl, ul et 77 sont declares avec l'attribut INTENT(IN), la 
procedure init ne peut done pas modifier leur valeur. II est interdit d'ecrire : rl = ul = ou 
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Exemple VII.2 

SUBROUTINE LECTURE (rr, uu, TT) 

IMPLICIT NONE 

REAL, INTENT (OUT) :: rr, uu, TT 

! Lecture des conditions initiales. 

OPEN (unit=12, f ile="condinit .dat" , status="old" , & 
form=" formatted" ) 

READ ( 12 , * ) 

READ (12, ' (llx,e8.3) ' ) rr 

WRITE (6,*) 'Density in g/cm-3: ', rr 
READ (12, ' (llx,e8.3) ' ) uu 

WRITE (6,*) 'Shock speed in cm/s: ', uu 
READ (12, ' (llx,f7.3) ' ) TT 

WRITE (6,*) 'T in eV: ', TT 
CLOSE (12) 

END SUBROUTINE LECTURE 

Dans l'exemple precedent, les arguments rr, uu et TT sont declares avec l'attribut intent(out), la 
procedure lecture doit done leur fournir une valeur. 

Remarque : cependant, l'experience montre que le comportement du compilateur depend aussi de 
son constructeur : certain compilateur ne signale aucune erreur, ni meme message d' attention, lors 
d'une mauvaise utilisation des variables a vocation INTENT(OUTJ. 

VII.A.2 Presence optionnelle des arguments 

L'attribut OPTIONAL permet de declarer certains arguments comme optionnels et leur presence 
eventuelle est testee a l'aide de la fonction intrinseque d'interrogation PRESENT. Cette precaution evite 
1' utilisation systematique de 1' argument, or que Ton sait pertinemment qu'il ne sera pas forcement 
toujours present. 

Syntaxe generate : 

type, optional [, attributs] :: var 
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Exemple VII.3 

SUBROUTINE LECTURE (rr, uu, TT, AtN) 
IMPLICIT NONE 

REAL, INTENT (OUT) :: rr, uu, TT 
REAL, OPTIONAL, INTENT (OUT) :: AtN 

! Lecture des conditions initiales. 

OPEN (unit=12, f ile="condinit .dat" , status="old" , & 
form=" formatted" ) 

READ ( 12 , * ) 

READ (12, ' (llx,e8.3) ' ) rr 

WRITE (6,*) 'Density in g/cm-3: ', rr 
READ (12, ' (llx,e8.3) ' ) uu 

WRITE (6,*) 'Shock speed in cm/s: ', uu 
READ (12, ' (llx,f7.3) ' ) TT 

WRITE (6,*) 'T in eV: ', TT 

IF ( PRESENT ( AtN ) ) THEN 

READ (12, ' (llx,i2) ' ) AtN 

WRITE (6,*) 'Atomic number: ', AtN 

END IF 
CLOSE (12) 

END SUBROUTINE LECTURE 

VII. A.3 Passage d' arguments par mots cles 

Dans un appel de procedure, il devient possible de reperer les arguments, non plus par leur 
position, mais aussi par le nom de l'argument correspondant dans l'interface de procedure. Cet appel 
par mot cle est tres pratique, surtout lorsque la liste des arguments contient des arguments optionnels. II 
est meme recommande d'utiliser les mot cles pour les arguments optionnels. Toutefois les arguments, 
qui ne sont encore reperes uniquement par leur position, doivent obligatoirement figurer en tete de liste, 
et seul le dernier pourra etre omis (s'il est optionnel). 

Exemple VII.4 

REAL SUBROUTINE EVALDIM (a, b, n, min, max, prec ) 

REAL, INTENT (IN), DIMENSION ( n ) : : a, b 

REAL, INTENT (OUT), OPTIONAL :: min, max, prec 

END SUBROUTINE EVALDIM 
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! Appels corrects 

CALL EVALDIM (x, y) 

CALL EVALDIM (x, y, min=eps ) 

CALL EVALDIM (min=eps, max=hug, a=x, b=y) 

CALL EVALDIM (min=eps, prec=mic, a=x, b=y) 

CALL EVALDIM (x, y, eps, hug) 

CALL EVALDIM (x, y, max=hug, prec=mic) 

VII.B INTERFACE DE PROCEDURE 

Une interface de procedure est constitute des informations necessaires permettant la 
communication entre plusieurs unites de programmes. Ces informations sont : 

- le type de la procedure : function ou subroutine ; 

- les arguments de la procedure (arguments formels), avec leur type et leurs attributs, 

- les proprietes du resultat dans le cas d'une fonction. 

Par defaut, on est dans le cas d'une interface implicite : les arguments affectifs de l'appel de la 
procedure sont definis dans 1' unite de programme appelant, et les arguments muets de la procedure sont 
definis dans la procedure. 

- Dans le cas, d'une procedure interne, la compilation de l'unite appelante et de la procedure se 
fait ensemble : il y a done un controle de coherence. C'est done une interface explicite, mais 
qui limite la visibilite de la procedure, depuis l'exterieur. 

- Dans le cas de procedures externes, la compilation de l'unite appelante et de la procedure se 
fait separement et les arguments sont passes par adresse. Alors, le controle de coherence entre 
arguments effectifs et arguments muet n'existe pas a la compilation. De nombreuses erreurs, 
lors de l'execution du programme, sont possibles et souvent difficiles a detectees. 

C'est pourquoi, on a recours a une interface explicite, a l'aide d'un bloc interface. Ainsi, le 
compilateur connait toutes les informations necessaires a l'interface entre l'unite appelante et la 
procedure. 

Syntaxe generate d'un un bloc interface : 
interface 

procedure nom_procedure (list_arg) 

bloc declaratif de list_arg 
end procedure nom_procedure 
end interface 
ou le bloc declaratif est une duplication de la partie declarative de la procedure. 
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VII.B.l Bloc interface dans l'unite appelante 

On ecrit directement le bloc interface au niveau du bloc declaratif de 1' unite de programme qui 
appelle la procedure externe concernee. 

Exemple VII.5 

PROGRAM CHOC 

IMPLICIT NONE 

REAL :: rrl, uul, TT1, rr2 , uu2 , TT2 
REAL, EXTERNAL : : FuncGas 

INTERFACE 

SUBROUTINE STRUCT (rl, ul, Tl, r2 , u2 , T2 , FuncGas, TG) 
REAL, INTENT (IN) :: rl, ul, Tl, r2 , u2 , T2 
REAL, EXTERNAL : : FuncGas 
CHARACTER ( LEN=2 ) , INTENT (IN) :: TG 

END SUBROUTINE STRUCT 
END INTERFACE 

CALL STRUCT (rrl, uul, TT1, rr2 , uu2 , TT2 , FuncGas, 
'GP' ) 

END PROGRAM CHOC 
SUBROUTINE STRUCT (rl, ul, Tl, r2 , u2 , T2 , FuncGas, TG) 
IMPLICIT NONE 



REAL, INTENT (IN) :: rl, ul, Tl, r2 , u2 , T2 
REAL, EXTERNAL : : FuncGas 
CHARACTER ( LEN=2 ) , INTENT (IN) :: TG 

INTEGER, PARAMETER :: nvar = 3, nstep = le4 
REAL, DIMENSION (nvar, 2) :: f 
INTEGER :: n, i, j, k, jm, km, count 
REAL :: rr, uu, TT, z, zO, h 
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Integre_choc : DO n = 1, nstep 

rr = f(l,l) * rl ; TT = f(3,l) * Tl 

IF ( rr>rho ( nmax ) .or. TT>TeV ( mmax ) ) THEN 

EXIT Integre_choc 
END IF 

CALL RUNGEKUTTA (z, h, f, 3, n, FuncGas ) 
END DO Integre_choc 

END SUBROUTINE STRUCT 

L' inconvenient, dans ce cas, est qu'il est necessaire de dupliquer le bloc interface dans toutes les 
unites de programme appelant la procedure concernee. Or l'interet d'une procedure externe, avec bloc 
interface, est justement d'etre utilisable, avec coherence, par toute autre unite de programme. On 
privilegiera alors l'ecriture du bloc interface au sein d'un module. 

VII.B.2 Bloc interface dans un module 

Pour ameliorer la fiabilite generale, on prefere inserer le meme bloc interface dans toutes les unites 
de programme faisant appel a la procedure en question. On utilise le module, qui via l'instruction use 
inserera le bloc interface dans l'unite appelante. 

A ce niveau, on a le choix entre deux voies : 

- ecrire un module dans le seul but de contenir le bloc interface ; 

- ecrire un module qui contient la procedure externe : c'est une excellente solution ! 

a) Module avec bloc interface 

On reprend l'exemple precedent, mais la declaration du bloc interface se fait grace a un module, 
nomme interhydro. 

Exemple VII.6 

MODULE INTERHYDRO 
INTERFACE 

SUBROUTINE STRUCT (rl, ul, Tl, r2 , u2 , T2 , & 

FuncGas, TG) 
REAL, INTENT (IN) :: rl, ul, Tl, r2 , u2 , T2 
REAL, EXTERNAL :: FuncGas 
CHARACTER ( LEN=2 ) , INTENT (IN) :: TG 
END SUBROUTINE STRUCT 
END INTERFACE 
END MODULE INTERHYDRO 
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PROGRAM CHOC 

USE INTERHYDRO ! Insertion de 1' interface 
IMPLICIT NONE 

REAL :: rrl, uul, TT1, rr2 , uu2 , TT2 
REAL, EXTERNAL : : FuncGas 

CALL STRUCT (rrl, uul, TT1, rr2 , uu2 , TT2 , FuncGas, 
'GP' ) 

END PROGRAM CHOC 

SUBROUTINE STRUCT (rl, ul, Tl, r2 , u2 , T2 , FuncGas, TG) 
IMPLICIT NONE 

REAL, INTENT (IN) :: rl, ul, Tl, r2 , u2 , T2 
REAL, EXTERNAL : : FuncGas 
CHARACTER ( LEN=2 ) , INTENT (IN) :: TG 

INTEGER, PARAMETER :: nvar = 3, nstep = le4 
REAL, DIMENSION (nvar, 2) :: f 
INTEGER :: n, i, j, k, jm, km, count 
REAL :: rr, uu, TT, z, zO, h 

Integre_choc : DO n = 1, nstep 

rr = f(l,l) * rl ; TT = f(3,l) * Tl 

IF ( rr>rho ( nmax ) .or. TT>TeV ( mmax ) ) THEN 

EXIT Integre_choc 
END IF 

CALL RUNGEKUTTA (z, h, f, 3, n, FuncGas) 
END DO Integre_choc 

END SUBROUTINE STRUCT 

b) Module contenant la procedure 

On reprend a nouveau l'exemple precedent, mais cette fois-ci le module sert a encapsuler la 
procedure, qui devient une procedure interne au module grace au mot cle contains. L'interface de 
procedure est connue vie l'instruction use nom_module. 
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Exemple VII.7 

MODULE INTERHYDRO 

CONTAINS 

SUBROUTINE STRUCT (rl, ul, Tl, r2 , u2 , T2 , FuncGas, TG) 
IMPLICIT NONE 

REAL, INTENT (IN) :: rl, ul, Tl, r2 , u2 , T2 
REAL, EXTERNAL : : FuncGas 
CHARACTER ( LEN=2 ) , INTENT (IN) :: TG 

INTEGER, PARAMETER :: nvar = 3, nstep = le4 
REAL, DIMENSION (nvar, 2) :: f 
INTEGER :: n, i, j, k, jm, km, count 
REAL :: rr, uu, TT, z, zO, h 

Integre_choc : DO n = 1, nstep 

rr = f(l,l) * rl ; TT = f(3,l) * Tl 

IF ( rr>rho ( nmax ) .or. TT>TeV ( mmax ) ) THEN 

EXIT Integre_choc 
END IF 

CALL RUNGEKUTTA (z, h, f, 3, n, FuncGas) 
END DO Integre_choc 

END SUBROUTINE STRUCT 
END MODULE INTERHYDRO 
PROGRAM CHOC 

USE INTERHYDRO ! acces directe a la procedure 
IMPLICIT NONE 

REAL :: rrl, uul, TT1, rr2 , uu2 , TT2 
REAL, EXTERNAL : : FuncGas 

CALL STRUCT (rrl, uul, TT1, rr2 , uu2 , TT2 , FuncGas, 
'GP' ) 

END PROGRAM CHOC 
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VII.B.3 Recapitulatif sur l'interface explicite 



Claire Michaut 



a) L'interface explicite 

Une interface de procedure est dite explicite dans les cinq cas suivants : 

- appel de procedures intrinseques : elles ont toujours des interfaces explicites : 

- appel de procedure interne ; 

- presence d'un bloc interface dans la procedure appelante, 

- acces au module, via l'instruction use, contenant le bloc interface de la procedure externe ; 

- acces au module, via l'instruction use, contenant la procedure externe ; 

b) L'interface explicite obligatoire 

Le long du cours, on a deja mentionne des cas ou l'interface doit etre explicite, sinon le 
compilateur ne peut pas savoir si le programme est coherent ou non. On rencontre dix situations de ce 
type : 

- utilisation d'une fonction a valeur tableau, 

- utilisation d'une fonction a valeur pointeur, 

- utilisation d'une fonction a valeur chaine de caracteres dont la longueur est determinee 
dynamiquement, 

- utilisation d'une procedure generique, 

- passage en argument d'un tableau a profil implicite, 

- argument formel avec l'attribut pointer ou target, 

- passage d'arguments a mots cle, 

- argument optionnel, 

- surcharge ou definition d'un operateur, 

- surcharge de l'operateur d'affectation. 



VII.C INTERFACE GENERIQUE 

Le bloc interface peut servir aussi a regrouper une famille de procedure, sous un meme nom. La 
procedure sera appelee dans une unite de programme par son nom de famille, ou nom generique. La 
procedure qui sera reellement utilisee lors de 1' appel est choisie pas le compilateur en fonction des 
instructions du bloc interface. Le bloc interface doit alors posseder un nom, qui sera le nom generique 
des procedures concernees. 
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Syntaxe generate : 

interface nom_generique 

procedure noml (list_argl) 

bloc declaratif de noml 
end procedure noml 
procedure nom2 (list_arg2) 

bloc declaratif de nom2 
end procedure nom2 

procedure nomn (list_argn) 

bloc declaratif de nomn 
end procedure nomn 
end interface 

Note : en Fortran 95, on peut aussi nommer l'instruction end interface : 
end interface nom_generique 

Exemple VII.8 

INTERFACE moyenne 

FUNCTION moy_int (tab_int) 

INTEGER, DIMENSION ( : ) :: tab_int 

END FUNCTION moy_int 

FUNCTION moy_real (tab_real) 

REAL, DIMENSION ( : ) :: tab_real 
END FUNCTION moy_real 

FUNCTION moy_comp (tab_comp) 
COMPLEX, DIMENSION ( : ) :: tab_comp 
END FUNCTION moy_comp 
END INTERFACE moyenne 

Dans cet exemple, la fonction moyenne accepte en argument un tableau de rangl soit d'entiers, 
soit de reels, soit de complexes. La moyenne sera effectivement calculee en faisant appel a la fonction 
correspondant au type de l'argument donne. 

VII.D SURCHARGE ET CREATION D 'OPERATEURS 

La notion de surcharge des operateurs predefinis par le langage est propre aux langages orientes 
objets. On entend par surcharge d'un operateur : un elargissement de son champ d' application. II faut, 
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bien entendu, definir les nouvelles operations entre les objets. La surcharge d'un operateur impose de 
respecter sa nature et les regies de priorite precedemment definie dans le langage. 

Habituellement, un operateur retourne une valeur, construite en fonction des expressions donnees. 
On emploie alors des procedures de type FUNCTION pour surcharger les operateurs. Toutefois, une 
exception est faite pour l'operateur d'affectation « = » qui ne retourne aucune valeur, il faut done 
utiliser une procedure de type SUBROUTINE, pour surcharger cet operateur. 

On definit la surcharge d'un operateur grace a un bloc INTERFACE. 

Remarque : en general, la surcharge est definie dans un module. 

Remarque : Certains operateurs ont deja fait l'objet d'une surcharge au sein du langage. 

VII.D.1 Interface OPERATOR 

Pour surcharger un operateur (autre que l'operateur d'affectation), on utilisera un bloc d'interface, 
en lui ajoutant l'instruction OPERATOR. Le bloc interface a la meme syntaxe et la meme fonction que 
dans le cas des interfaces explicites. C'est le compilateur qui choisit quelle procedure sera effectivement 
utilisee lors de tel ou tel appel. 

On indique l'operateur que Ton veut surcharger, entre parentheses, a la suite du mot cle 
OPERATOR : 

- si, c'est un operateur existant, on utilise directement son signe ; 

- si, on cree un nouvel operateur, on utilise un identifiant, place entre les caracteres « . ». 
Syntaxe generate : 

interface operator (*) (ou autre signe) 

function nom_fun (list_arg) 

bloc declaratif avec intent (in) 
end function nom_fun 

end interface 

ou : 

interface operator (.nom_op.) 
function nom_fun ( list_arg) 

bloc declaratif avec intent (in) 
end function nom_fun 

end interface 

Remarque : Les arguments de la fonction, associee a un operateur pour sa sur-definition, doivent 
avoir l'attribut intent (in). 
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Exemple VII.9 

MODULE TYPMAT 

TYPE MAT_INT 

INTEGER : : lin, row 

INTEGER :: DIMENSION (:,:) , POINTER :: ptr_mat 
EN TYPE MATRICE 
END MODULE TYPMAT 



INTERFACE OPERATOR (*) 
FUNCTION MULT (a, b) 
USE TYPMAT 

TYPE (mat_int), INTENT (IN) :: a, b 
TYPE (mat_int) :: mult 
END FUNCTION MULT 
END INTERFACE 



TYPE (mat_int) FUNCTION MULT (a, b) 
USE TYPMAT 

TYPE (mat_int) , INTENT (IN) :: a, b 
INTEGER, TARGET, PRIVATE :: i, j 
INTEGER : : ok 

i = a%lin ; j = b%row 
mult%lin = i ; mult%row = j 

ALLOCATE (mult%ptr_mat(i, j ) , stat=ok) 
IF (ok > 0) THEN 

PRINT*, "Erreur d' allocation de mult" 

STOP 120 
END IF 

mult%ptr_mat = MATMUL ( a%ptr_mat , b%ptr_mat) 
END FUNCTION MULT 

VII.D.2 Interface ASSIGNMENT 

Comme nous l'avons mentionne precedemment, l'operateur d' affectation se comporte 
differemment des autres operateurs, car il ne retourne pas une nouvelle valeur. Done pour surcharger cet 
operateur, on utilise 1' instruction ASSIGNMENT (=). La syntaxe est identique au bloc INTERFACE 
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OPERATOR, se reporter done a l'exemple precedent. Par contre, on utilisera une procedure de type 
subroutine pour sur definir 1' affectation. 

Syntaxe generate : 

interface assignment (=) 

subroutine nom_fun (list_arg) 

bloc declaratif avec intent (out ) 

puis intent (in) 
end subroutine nom_fun 

end interface 

Remarque : Les arguments de la subroutine, associee a l'operateur « = » pour sa sur-definition, 
doivent avoir respectivement : 

- l'attribut intent (out) ou intent (inout), pour le premier argument, qui est l'operande de 
gauche ; 

- l'attribut intent ( in), pour le second argument, qui est l'operande de droite. 
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VIII LES ENTREES / SORTIES 

Le langage Fortran dispose d'un ensemble de possibilites en matiere d'entree / sortie (E/S) 
particulierement riches. 

Les nouvelles caracteristiques significatives apportees par la norme Fortran 90 sont: 

- les E/S sans deplacement (nom-advancing I/O), 

- les listes nommees (NAMELIST), 

- quelques nouveaux specificateurs des instructions OPEN et INQUIRE. 

- quelques nouveaux descripteurs d'edition et une generalisation du descripteur d'edition G, 

VIII.A INSTRUCTIONS GENERALES 

Les instructions generates d'entree / sortie sont les instructions de lecture et d'ecriture : READ, 
WRITE et PRINT. On distingue les entrees / sortie a acces directes des entrees / sortie a sequentielles. 

- L'acces sequentiel permet de traiter les enregistrements dans l'ordre dans lequel ils 
apparaissent ; ils sont done reperes par la position du pointeur. 

- Tandis que l'acces direct gere des enregistrements de taille unique dans un ordre quelconque ; 
ils sont done reperes par un numero d'enregistrement (REC). 

Syntaxe generate acces sequentiel : 
read ( [unit=]unit, [ fmt= ] format [, iostat] [, err] 

[, end] [, advance] [, size] [, eor]) list_var 
write ([unit=] unit, [ fmt= ] format [, iostat] [, err] 

[, advance]) list_var 
print fmt, list_var 

Syntaxe generate acces direct: 
read ( [unit=]unit, [fmt=] format [, rec ] [, iostat] 

[, err] [, end]) list_var 
write ( [unit=]unit, [fmt=] format [, rec] [, iostat] 
[ , err] ) list_var 

Les parametres sont les suivants : 

- unit designe 1' unite dans laquelle on veut lire ou ecrire ; 

- fmt est l'indicateur de format ; 

- rec designe le numero de l'enregistrement ; 

- iostat indicateur d'erreurs ; 

- err renvoie a l'etiquette d'instruction en cas d'erreurs ; 

- end renvoie a l'etiquette d'instruction lorsque la fin de fichier est rencontree ; 
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- advance precise si le pointeur avance ou pas a la fin de l'enregistrement ; 

- size permet de recuperer la longueur d'enregistrement ; 

- eor renvoie a l'etiquette d'instruction lorsque la fin d'enregistrement est rencontree. 

VIII. A. 1 Description des parametres de lecture / ecriture 

a) l'unite 

Lorsqu'on veut acceder a un fichier, on donne au parametre unit la valeur entiere qui repere le 
fichier. Le numero d'identification du fichier est defini a son ouverture (voir l'instruction OPEN ci- 
apres). 

Lorsqu'on veut acceder a une variable du type chaine de caracteres, on affecte au parametre unit le 
nom de la variable. 

Exemple VIII. 1 

INTEGER :: n, k 

CHARACTER (LEN=10), PRIVATE :: name = ' graphX . out ' 
INTEGER, PRIVATE : : iname = 

DO k = 1, n 

iname = iname + 1 

WRITE (name(6:6) , ' (il) ' ) iname 
EN DO 

b) le format 

Pour ecrire un format d'enregistrement, il faut connaitre quelques regies : on precise le type 
d'enregistrement, le nombre total de caracteres, le nombre de chiffres significatifs pour les nombres 
decimaux..., puis on peut affiner le format avec des espaces. . . 

Les descripteurs de types : 

- a : pour une chaine de caracteres ; 

- 1 : pour une valeur logique ; 

- i : pour un nombre entier ; 

- f : pour un nombre decimal ; 

- e : pour un reel en notation scientifique ; 

Ensuite, on donne le nombre de caracteres directement : a3, par exemple. 

Lors de l'utilisation de formats identiques, les uns a la suite des autres, on peut multiplier les 
formats directement avec un nombre de fois : 4a3, par exemple, pour ecrire 4 fois une chaine de 3 
caracteres. Pour des format plus complexes, on se sert de parentheses et on opere une distribution 
comme pour la multiplication. Par exemple : 4(a3, i2) 
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On peut sauter des caracteres, avec l'identifiant x : 4(a3, lx, i2). Ainsi, on saute 1 caractere entre la 
chaine de 3 caracteres et l'entier a 2 chiffres ; le tout est repete 4 fois. En ecriture, cela revient a laisser 
un espace. 

On peut se positionner dans la memoire tampon, avec le descripteur t. Pour cela, on indique aussi 
Femplacement sur lequel on veut pointer. Par exemple : t20, veut dire se placer sur le 20 ieme caractere de 
la ligne. 

Le caractere « / » indique d'aller a une nouvelle ligne. 

On peut aussi appliquer un facteur d'echelle, sur les nombres reels, a l'aide de l'identifiant p. Par 
exemple, 3p devant les autres descripteurs, multiplie par 10 3 la mantisse et diminue l'exposant de 3. 

Ces descripteurs de format sont les plus couramment utilises, mais la liste des descripteurs n'est 
pas, ici, exhaustive, ni meme leurs effets. Pour plus de details, il faut se referer a un manuel complet du 
langage Fortran. 

Exemple VIII.2 

FORMAT ("vail = ", lp, f7.4, lx, "val2 = ", i3 ) 

FORMAT ("Les resultats sont : ", 3(el2.4, 2x, f8.3)) 

FORMAT (tl2, 3(el2.4, 2x, f8.3), 3x, i20) 

c) le numero de l'enregistrement 

Ce parametre (rec) ne sert que dans les fichiers a acces direct et repere le numero de 
l'enregistrement. C'est done une valeur entiere. 

d) rindicateur d'erreurs 

On utilise le parametre iostat, pour gerer les erreurs d'enregistrement. Ce parametre peut prendre 
plusieurs valeurs, de type entier : 

- il est nul, lorsqu'il n'y a aucune erreur ; 

- il est negatif lorsque l'erreur est generee par un fin : 

- de fichier {end), 

- ou d'enregistrement (eor), dans les acces sequentiels formates, sans 
avancement du pointeur (advance = 'no') ; 

- il est positif dans les autres cas d'erreurs. 

e) err renvoie a l'etiquette destruction en cas d'erreurs 

Les trois parametres err, end et eor recoivent chacun une etiquette, qui renvoie a l'instruction 
correspondante, en fonction des besoins. Ainsi, le programmateur peut decider de la suite du 
programme, lorsque des erreurs surgissent dans la lecture ou l'ecriture d'enregistrement, lorsqu'une fin 
de fichier, ou bien une fin d'enregistrement, sont respectivement rencontrees. 
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Exemple VIII.3 : etiquettes d'instruction 

CHARACTER (LEN=4) :: nom 
DO 

READ (15, 'i4', err=105, end=106) nom 
END DO 

105 PRINT*, 'erreur a la lecture' ; STOP 222 

106 PRINT*, 'on a fini !' 

Note : Le parametre iostat permet aussi d'orienter le programme si une erreur apparait, ou si une 
fin de fichier est detectee. L'interet est de permettre une programmation plus structured que celle de 
l'exemple precedent. 

Exemple VIII.4 

CHARACTER ( LEN=4 ) :: nom 

INTEGER : : dem 

DO 

READ (15, ' i4', iostat=dem) nom 

IF (dem > 0) THEN 

PRINT*, 'erreur a la lecture' ; STOP 222 

ELSE 

PRINT*, 'on a fini !' ; EXIT 
END IF 
END DO 

f) avancement ou non du pointeur 

Par defaut, lors de l'execution d'une instruction d'entree / sortie, le pointeur d'un fichier passe a 
Fenregistrement suivant. Toutefois, on peut preciser que le pointeur reste sur le dernier enregistrement, 
avec le parametre advance. 

- advance = 'no ' specifie que le pointeur reste sur place ; 

- advance = 'yes' specifie que le pointeur avance a Fenregistrement suivant (c'est le defaut). 
Dans le cas d'enregistrement statique (advance = 'no'), le format doit toujours etre fixe. 



73 



Ecole doctorale d'Astronomie et d'Astrophysique d'tle-de-France - Cours de Fortran 90 1 95 - Claire Michaut 

g) longueur d'enregistrement 

Dans le cas ou une fin d'enregistrement est detectee, le parametre size permet de recuperer la 
longueur de l'enregistrement. Ce parametre ne s' applique qu' a une instruction de lecture statique 
(advance = 'no'), et done formatees. 

READ (12, fmt=405, advance = 'no', size = 1, eor=106) var 

106 PRINT* , 'longueur ', 1 



VIII.B OUVERTURE DE FICHIERS 

L'instruction OPEN permet d'ouvrir un fichier, la syntaxe generale est la suivante : 

open (unit, file, status, action, access, iostat, err, 
position, form, reel, blank, pad, delim) 

- unit : numero d'identification de l'unite logique que Ton manipule ; 

- file : nom du fichier que Ton ouvre ; 

- status :statut du fichier (nouveau, ancien. . .) ; 

- action : precise ce qu'on veut faire de ce fichier (lecture, ecriture) ; 

- access : precise l'acces direct ou sequentiel a un fichier ; 

- iostat : sert a mentionner s'il y a des erreurs a l'execution ; 

- err : sert a donner une instruction en cas d'erreurs ; 

- position : indique la position du pointeur du fichier ; 

- form : precise si un fichier est formate ou non ; 

- reel : precise la longueur d'un enregistrement ; 

- blank : precise la gestion des espaces lors de la lecture ; 

- pad : sert pour completer ou non avec des blancs les enregistrements ; 

- delim : sert a delimiter des chaines de caracteres. 



VIII.B. 1 Description des parametres d'ouverture 

a) l'unite logique 

On repere un fichier par un numero, en affectant une valeur au parametre unit. Tant que le fichier 
reste ouvert, il est identifie par ce numero dans les instructions d'entree / sortie du reste programme. 

Exemple VIII.5 

OPEN (unit=2 0) 

READ (20, fmt=12) var 

b) nom du fichier 

On affecte, au parametre file, une chaine de caracteres qui est le nom du fichier a ouvrir. 
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Exemple VIII.6 

OPEN (unit=20, f ile= ' f ich_res ' ) 

c) status 

Ce parametre peut recevoir une de ces cinq affectations : new, old, unknown, replace, scratch. 

- new : si le fichier doit etre cree. 

- old : si le fichier doit deja exister. 

- unknown : statut inconnu qui depend de l'environnement. 

- replace : si le fichier existe deja lors de l'ouverture, il sera detruit et dans tous les cas un 
nouveau fichier sera cree. 

- scratch : indique qu'un fichier anonyme sera cree et connecte a l'unite specifiee via le 
parametre unit. La duree de vie de ce fichier sera soit le temps d'execution du programme, soit 
le temps s'ecoulant jusqu'a la fermeture de cette unite. 

d) vocation du fichier 

Le parametre action peut recevoir les attributs read, write, readwrite. 

- read : toute tentative d'ecriture est interdite. 

- write : toute tentative de lecture est interdite. 

- readwrite : les operations de lecture et d'ecriture sont autorisees. 

e) acces direct ou sequentiel 

Pour preciser le type d'acces a un fichier, on affecte la parametre access des attributs sequentiel ou 
direct. L'acces sequentiel etant Faeces par defaut. 

Exemple VIII.7 

OPEN (unit=20, f ile= ' f ich_res ' , action^ 'write ', & 
acces= 'direct ' ) 

f) detection d'erreurs 

On se sert du parametre iostat pour detecter les erreurs lors de l'ouverture : 

- iostat = : il n'y a pas d'erreurs : 

- iostat > : des erreurs sont survenues. 

g) instruction en cas d'erreurs 

On affecte une etiquette au parametre err. En cas d'erreur detectees, le programme execute les 
instructions contenues a partir de l'etiquette. (Voir Exemple VIII.3 : etiquettes d'instruction.) 

h) position du pointeur du fichier 

Le parametre position peut recevoir les attributs rewind, append, asis. 

- rewind : le pointeur du fichier sera positionne au debut du fichier. 
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- asis : permet de conserver la position du pointeur de fichier. II est utile lorsqu'on veut modifier 
certaines caracteristiques du fichier, tout en restant positionne au meme endroit. 

- append : le pointeur du fichier sera positionne a la fin du fichier. 

i) type formate ou non 

On specifie le type des enregistrements, contenus dans un fichier, a l'aide du parametre form, qui 
recoit ainsi les attributs formatted ou unformatted. Par defaut, les enregistrements sont formates. 



j) longueur d'un enregistrement 

La parametre reel specifie la longueur d'un enregistrement pour l'acces direct, sinon precise la 
longueur maximale d'un enregistrement pour l'acces sequentiel. 



k) blank 

Ce parametre ne s'utilise que dans le cas d' enregistrements formates. II peut prendre les attributs 
null, ou zero. 

- null : on ne tient pas compte des espaces, comme si il n'existaient pas. 

- zero : on remplace les espaces par des 0. 

1) pad 

Ce parametre recoit les attributs yes, ou no. 

- yes : un enregistrement lu avec format est complete par des blancs, si necessaire. 

- no : pas de complement avec des blancs. 

m) delimitation des chaines de caracteres 

Le parametre delim peut recevoir les attributs apostrophe, quote, none. 

- apostrophe : indique que l'apostrophe sera utilise pour delimiter les chaines de caracteres lors 
d'une ecriture de type namelist. 

- quote : indique que l'apostrophe sera utilisee pour delimiter les chaines de caracteres lors d'une 
ecriture de type namelist. 

- none : indique qu'aucun delimiteur ne sera utilise (valeur par defaut) 



VIII.C L'INSTRUCTION INQUIRE 

C'est un instruction d'interrogation sur les parametres d'un fichiers. Elle renvoie done une valeur 
logique « .TRUE. » ou « .FALSE ». Elle possede des parametres tres proches de ceux de l'instruction 
open : position, delim, action, pad ont la meme signification et les memes valeurs que pour l'instruction 
open. 
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Syntaxe generate : 

inquire (unit, file, iostat, err, 

opened, exist, number, named, name, 

read, write, readwrite, action, access, direct, 

position, form, reel, nextrec, blank, pad, delim) 

a) read, write, readwrite 

Ces parametres recoivent chacun les attributs yes, no ou unknown. ; ils indiquent respectivement 
si la lecture, l'ecriture et la lecture et l'ecriture, dans un fichier, sont autorisees ou non, ou 
indeterminees. 

VIII.C.2 Par liste de sortie 

Une variante de l'instruction INQUIRE permet de se renseigner sur la longueur d'une liste de sortie 
lorsque l'ecriture d'enregistrements non formates est souhaitee. Cette variante est connue sous le nom 
"INQUIRE par liste de sortie", elle a la syntaxe suivante : 

inquire (iolength = longueur) out_list 

ou longueur est une variable scalaire de type entier par defaut, utilisee pour determiner la longueur de la 
liste de sortie non formatee outjtist. 

Note : cette longueur peut etre utilisee pour determiner la valeur a fournir au specificateur reel lors 
d'une prochaine instruction open. 

Exemple VIII.8 

INTEGER : : longueur 

INQUIRE ( IOLENGTH=longueur ) nom, adresse, tel 

OPEN (1, FILE= ' repertoire ' , RECL=longueur , & 
FORM= ' unformatted ' ) 

WRITE(l) nom, adresse, tel 

VIII.D L'INSTRUCTION NAMELIST 

L'instruction namelis t permet de lier un ensemble de variables dans un seul groupe, afin de 
faciliter les operations d'entree /sortie. 
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VIII.D.1 Declaration 



La declaration d'un tel groupe se fait par une instruction de specification namelist avant toute 
instruction executable, dont la forme est : 

namelist /nom_namelist/ liste_var 

Remarques : 

- II est possible de completer Enumeration de la liste dans la meme unite de portee. 

- Si une variable a l'attribut PUBLIC, aucune variable dans la liste nommee ne peut avoir 
l'attribut PRIVATE. 

- Les champs doivent etre des variables definies a la compilation, done pas de parametres de 
tableaux ou de chaines, d'objets automatiques ou pointers. . . 

VIII.D.2 Usage 

Dans une instruction READ ou WRITE, la liste nommee sera designee par le specificateur nml . 
Exemple VIII.9 

INTEGER, DIMENSION(5) :: age = (/6, 7, 8, 9, 10/) 
CHARACTER ( LEN=3 ) :: classe ( 5 )=( / ' cp ' , ' eel ' , ' ce2 ' , 

' cml ' , ' cm2 ' / ) 
CHARACTER (LEN=15) :: maitre(5) = " 
NAMELIST /ecole/ age, classe, maitre 

• • • 

READ ( * , nml=ecole) 

! lit des donnees sous la forme suivante: 
&ECOLE maitre(4) = 'Mme Grosso' / 
&ECOLE maitre(2) = 'M. Chaduiron' / 

WRITE(*, nml= ecole) 

! genere les lignes suivantes : 

&ECOLE AGE = 6, 7, 8, 9, 10, CLASSE = cpcelce2cmlcm2 
MAITRE = M. Chaduiron Mme Grosso 

/ 
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IX ANNEXE 



IX.A QUELQUES FONCTIONS INTRINSEQUES 
IX. A. 1 Manipulation de bits 

Les variables i, j, n, m, utilisees sont toutes des valeurs entieres. On rappelle que les bits d'un 
nombre sont numerates de gauche a droite et comment avec le bit 0. 

iand (i, j) 
ieor (i, j) 
ior (i, j) 

Ces trois fonctions renvoient respectivement un entier du meme nombre de bits que i, forme par le 
resultat des trois operations logiques : ET, OU exclusif, OU inclusif, effectuees bit a bit entre les entiers 
ietj. 

ishft (i, shift [, size]) 
ishftc (i, shift [, size]) 

Ces fonctions precedent a un decalage vers la gauche, ou vers la droite, du nombre de bits indique 
par shift. Lorsque shift est positif, on decale vers la gauche, lorsque shift est negatif, on decale vers la 
droite. La fonction ishft remplie les bits manquants par des 0, tandis que la fonction ishftc opere une 
permutation circulaire. Elles renvoient toutes les deux un entier du meme nombre de bits que l'entier i. 

Exemple IX. 1 

ishft(120,2) donne 480 

ishftc ( 134678, 8) donne 34477568 

not (i) 

Cette fonction fait le complement a 1 (ou complement logique) des bits de l'entier i ; elle retourne 
un entier du meme nombre de bits que i. 

btest (i, n) 

Cette fonction renvoie « .TRUE. » si le nieme bit du nombre entier i est a 1, sinon « .FALSE. » si 
il est a 0. 
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ibset (i, n) 
ibclr (i, n) 

Ces deux fonctions retournent un entier identique a i, avec respectivement le nieme bit a 1 et le 
nieme bit a 0. 

ibits (i, n, m) 

Cette fonction extrait m bits de Tender i, a partir du nieme bit ; puis les decale vers la droite pour 
former un entier du meme type que i (les autres bits de gauche sont mis a 0). 

mvbits (i, n, m, j, nj ) 

mvbits est un sous-programme, dons on l'appelle avec un call), il est construit sur le meme 
principe que la fonction ibits. II extrait m bits de Tender i, a partir du nieme bit ; puis les place a la 
« njieme » position de Tender j. 

IX.A.2 Precision et codage numerique 
digit (x) 

Cette fonction admet un argument x entier ou reel et renvoie un entier egal au nombre de chiffres 
significatifs de x: 

- le nombre de chiffres binaires de la mantisse si x est reel 

- le nombre de chiffres binaires de stockage (hors bit de signe) si x est un entier. 

epsilon (x) 

Cette fonction renvoie un reel avec meme parametre de sous-type que x qui est presque 
negligeable compare a 1 . 

huge (x) 

Cette fonction admet un parametre x entier ou reel et renvoie la plus grande valeur representable 
dans le sous-type de x (limite d'overflow). 

tiny (x) 

Cette fonction admet un parametre x reel et renvoie la plus petite valeur positive non nulle 
representable dans le sous-type de x (limite d'underflow). 

maxexponent (x) 

Cette fonction admet un parametre x reel et renvoie l'entier representant le plus grand exposant 
possible dans le sous-type de x. 
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minexponent (x) 

Cette fonction admet un parametre x reel et renvoie l'entier representant le plus petit exposant 
possible dans le sous-type de x. 

exponent (x) 

Cette fonction renvoie un entier egal a l'exposant de la representation numerique de x (0 si x est 
egal a 0) 

fraction (x) 

Cette fonction renvoie un reel qui est egal a la mantisse de la representation numerique de x. 
nearest (x, s) 

Cette fonction renvoie un reel ayant le meme parametre de sous-type que x, egal au nombre 
exactement representable en machine le plus proche de x dans la direction indiquee par l'argument s. 

spacing (x) 

Cette fonction renvoie un reel ayant le meme parametre de sous-type que x, egal a l'espacement 
absolu des nombres de meme parametre de sous-type que x dans le voisinage de x. 

range (x) et precision (x) 

Ces fonctions sont decrites au chapitre II.C.4c) range et precision. 

IX.A.3 Fonctions numeriques elementaires 

a) Avec conversion 
ceiling (a) 

Cette fonction renvoie le plus petit entier du type par defaut qui est superieur ou egal a a. 
floor (a) 

Cette fonction renvoie le plus grand entier du type par defaut qui est inferieur ou egal a a. 

b) Sans conversion 
modulo (a,p) 

Cette fonction renvoie a modulo p pour a et p soit tous les deux reels soit tous les deux entiers. 
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IX.A.4 Fonctions elementaires de type caractere 

a) Conversion entier-caractere 
achar (i) 

Cette fonction renvoie le caractere dont le code ASCII est specifie par l'argument i. 
iachar (c) 

Cette fonction renvoie le code ASCII correspondant au caractere c. 
ichar (c) 

Cette fonction renvoie un entier qui represente le numero du caractere c dans la table a laquelle 
appartient c (ASCII/EBCDIC) 

b) Manipulation de chaine de caracteres 
adjustl (string) 

Cette fonction cale a gauche la chaine de caractere STRING. 

adjustr (string) 

Cette fonction cale a droite la chaine de caractere STRING. 

len_trim 

Cette fonction renvoie un entier qui represente la longueur de chaine de caractere STRING (sans 
compter les caracteres blancs a la fin). 

scan (string, set [,back]) 

Cette fonction renvoie un entier qui represente le numero du premier caractere de string ou 
apparait l'un des caracteres de set ou sinon. back est optionnel et permet de trouver la derniere 
occurrence. 



verify (string, set [,back]): 

Cette fonction renvoie un entier qui represente le numero du premier caractere de string ne 
figurant pas dans set ou sinon. back est optionnel, s'il est egal a « .TRUE. », la recherche se fait depuis 
la fin de string. 
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c) Fonctions de comparaison lexicale 
lge (stringa, stringb) 
lie (stringa, stringb) 

Ces fonctions renvoient la valeur « .TRUE. », respectivement si stringa est apres ou avant stringb 
(ou au meme niveau) dans la table ASCII, sinon « .FALSE. ». 

lgt (stringa, stringb) 
lit (stringa, stringb) 

Ces fonctions renvoient la valeur « .TRUE. », respectivement si stringa est strictement apres ou 
avant stringb dans la table ASCII, sinon « .FALSE. ». 

Remarque : en cas d'inegalite de longueur, la chaine la plus courte sera completee a blanc sur sa 
droite. 

repeat (string, n) 

Cette fonction permet la formation d'une chaine de caracteres par concatenation de n fois la chaine 
string. 

trim( string) 

Cette fonction renvoie la chaine de caracteres string sans les blancs eventuels situes a sa fin. 

IX.B FONCTIONS DE GENERATION DE NOMBRES ALEATOIRES 
random_n umber (harvest) 

Ce sous-programme genere dans harvest un nombre (ou des nombres) pseudo aleatoire, dans 
l'intervalle [0,1 [ ; harvest doit etre de type reel et peut etre un tableau. 

random_seed ([size] [, put] [, get]) 

Ce sous-programme interroge ou modifie le generateur de nombres aleatoires. II n'admet qu'un 
seul des trois arguments a la fois. 

- size : type entier d'attribut OUT ; le sous-programme retourne alors la taille du tableau 
d'entiers utilises comme germe pour la generation des nombres aleatoires. 

- put : type entier d'attribut IN ; le sous-programme utilise alors ce tableau d'entiers et de rang 1 
comme germe pour la generation des nombres aleatoires. 

- get : type entier d'attribut OUT ; le sous-programme retourne alors le tableau d'entiers et de 
rang 1 utilise comme germe pour la generation des nombres aleatoires. 
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IX.C HORLOGE EN TEMPS REEL 

date_and_time (date, time, zone, values) 

Ce sous-programme retourne dans les variables caracteres date et time la date et l'heure en temps 
d'horloge murale. L'ecart par rapport au temps universel est fourni optionnellement par zone, values est 
un vecteur d'entiers recuperant les informations precedentes sous forme d'entiers. 

system_clock (count, count_rate, count_max) 

Ce sous-programme retourne dans des variables entieres la valeur du compteur de periode 
d'horloge : count, le nombre de periodes par seconde : count_rate et la valeur maximale du compteur de 
periode : countjmax. On peut, par exemple, evaluer le temps CPU consomme par une portion de 
programme. 
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ix.d Sites interessants 

FORTRAN 90 - procedures predefinies : http://www.univ-lillel.fr/~eudil/pyrz/home.htm 
Le Programmeur : Fortran 90 : http://www.stcommunications.ca/~mboisso/Fortran 90 
Bibliotheque Fortran-90 orientee Objet : http://www.math.u-psud.fr/~laminie/F90_lib 
Cours IDRIS Fortran90 / 95 : http://www.idris.fr/data/cours/lang/f90 
http://consult.cern.ch/cnl/215/node28.html 

Fortran 90/ 95 - Le site Francophone : http://www2.cnam.fr/~lignelet/ftn.htmlx 
http://cch.loria.fr/documentation/documents/F95 
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Programmer en Fortran 90 - Guide Complet 
Auteur : Claude Delannoy 
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